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Localized wave solutions of the scalar homogeneous wave 
equation and their optical implementation 

Kaido Reivelt and Peeter Saari 
Institute of Physics University of Tartu 

Abstract 

In recent years the topic of localized wave solutions of the homogeneous scalar wave equation, 
i.e., the wave fields that propagate without any appreciable spread or drop in intensity, has been 
discussed in many aspects in numerous publications. In this review the main results of this rather 
disperse theoretical material are presented in a single mathematical representation - the Fourier 
decomposition by means of angular spectrum of plane waves. This unified description is shown to 
lead to a transparent physical understanding of the phenomenon as such and yield the means of 
optical generation of such wave fields. 
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I. INTRODUCTION 



The birth of the diffraction theory of hght dates back to the works of Francis Maria 
Grimaldi (1618 - 1663), Robert Hook (1635 - 1703), Christiaan Huygens (1629 - 1695) and 
Thomas Young (1773 - 1829) and was mathematically formulated by Augustin Jean Fresnel 
(1788 - 1827). Over the two centuries it has been considered as a very successful theory - 
it indeed very precisely describes the propagation of light in linear media. The foundation 
of the diffraction theory is the principle of Huygens which states that (i) all the points of 
a wavefront act as the sources of secondary wavelets and (ii) the field at all the subsequent 
points is determined by the superposition of those wavelets. 



The topic of free-space propagation o 



when James Neill Brittingham claimed 



1 vv 



' wave fields has attracted a renewed interest in 1983 
that he discovered a family of three-dimensional, 
nondispersive, source-free, free-space, classical electromagnetic pulses which propagate in a 
straight line in free space at light velocity (in this work he also introduced the term focus 
wave mode (FWM) for those wave fields). Now, the very idea of secondary spherical sources 
in the classical diffraction theory implies that any optical wave field suffers from lateral 
and longitudinal spread in the course of propagation in free space, whereby the diffraction 
angle of the spread is the larger the narrower is the field radius. In the view of this general 
principle the Brittingham's statement is an astounding one and he quite rightly used the 
formula "to convince the scientific community" in its arguments. However, the original 
focus wave mode was indeed the solutions of the Maxwell's equations and the scientific 
community had to resolve this apparent contradiction. As to give the reader an idea of the 
initial problem the theoreticians had to tackle with, we reproduce here the original definition 
of Brittingham which he deduced by " a very extensive heuristical fit of various differential 
equation solutions": given the Maxwell equations (SI) 

V ■ D = 
V-B = 0, 

where E, D, H, B and t are electric field, electric flux density, magnetic field, magnetic 
induction and time variable, respectively and using cylindrical coordinate system (p, ip, z—ct) 



the mathematical formulation of the original FWM reads as 

Hp{p,<f,z,t) = ^'3 + ^; 

(p, if, z, = *4 + % 
H,{p,ip,z,t)^^5 + %, 

where the functions {q — 1, 2. ..5) are written as 

- Ag (p, z - ct) Gi (p, z - ct) G2 {z - ct) G3 {z - cit) (0) 
for g = 1 and 4, and 

= Aq {p, Z - ct) Gi {p, Z - ct) G2 {z - Ct) G3 {z - Cit) $ (0) 



for g = 2, 3, 5. In those equations 



Gi — exp 



PL 

4F 



G2 — exp [—iki {z — ct)] 
G3 = exp [ik2 {z — cit)] 



and 

F^ig{z-ct)+C . 

The remaining definitions read as 

D^'^ r(n + l)c^p"-i --"+1 



A, 



pn+2 



cgp 

^pn+3 



+ 



{kic — k2c) p" ^ 

P'n+2 



D 



TE 



n{n + l)cgp" {3n + A)cgp 



,n+l 



pn+2 



n-1 



+ 



n {k\c — /C2C1) p' 

pn+l 

2cgp''+^ 2 (/cic - /C2C1) p"+^ 



IQpn+i 
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A:, 



D 



TE 



n{n + l) gp^ ^ (3n + 4) gp 



n+l 



+ 



pn+2 

n {-ki + k2) p^-^ 

pn+l 

2c,p«+3 2 (-fci + k2) 



^4 



TE 



[n + 1) 5fp^ 

pn+2 



n-1 



+ 



,n+l 



(-A;i + A;2)p' 



,n-l 



4. = -iD 



TE 



P 



n+2 



(n + 1) p*^ 



4^n+3 pn+2 

where D^^, (7, ^, ki and /c2 are constants. The $ functions are defined as 

sin {rup) 
cos (nc/?) 

n cos (nv?) 
— n sin [nip) 




and the supplemental conditions read 



where 



2gk2di — 1 
9 



0^2 



1 - 



One has to agree, that the physical idea is very much hidden behind this mathematical 
formulation. 

Brittingham claimed, that this mathematical formulation (i) satisfy the homogeneous 
Maxwell's equations, (ii) is continuous and nonsingular, (iii) has a three-dimensional pulse 
structure, (iv) is nondispersive for all time, (v) move at light velocity in straight lines, and 
(vi) carry finite electromagnetic energy. Thus, the formulas above give a mathematical 
formulation of a free-space wave field that can be described as a "light bullet" and, though 



the proof of the last claim was shown to be faulty by Wu a nd K ing lal, the whole idea was 

n n 

very intricate and rose a considerable scientific interest |3l|- |149| . 

The theoretical work of following years could be divided into the following topics (see 
also Ref. l35| for an overview): 



In the following publications 0, 0] the original vector field was reduced to its scalar 
counterpart and the dominant part of the research work that followed has been formulated 
in terms of solutions to homogeneous scalar wave equation. 

The close connection between the FWM's and the solutions of the paraxial wave equation 
and Schrodinger's equation (which both allow localized solutions) has been established 13, 

nn 

|5|, |6|| - it has been shown that in terms of the variables z + ct and z — ct, if the solution of 
the scalar wave equation is given by the anzatz exp [/5 {z + ct)] F (x, y,z — ct), the problem 
can be reduced to one of those equations. 



The infinite energy content of the original FWM's has been addressed in several pub 



lications (see Refs. 

of all, Sezginer 0] and Wu and Lehmann 




I45 . 4fi[ 0| and references therein). First 
^ proved that any finite energy solution of the 



wave equation irreversibly leads to dispersion and to spread of the energy. Then Ziolkowski 
pointed out, that the superpositions of the infinite energy FWM's could result in fi- 
nite energy solutions and in following publications a number of finite energy solutions to 
the scalar wave equation and Maxwell equations were deduced - "electromagnetic directed- 
energy pulse trains" (EDEPT) 0, |3| , " acoustic directed-energy pulse trains" (ADEPT) 



41 1 , splash pulses 0|, modified power spectrum (MPS) pulses, electromagnetic missiles 



(M 

Q 



1^1 



2a |22|, various super- and subluminal pulses [l^ etc. In correspondence with this 
broader class of localized waves (LW) have generally extended but finite r ang es of localiza- 
tions. Also, several alternate infinite energy LW's (Bessel-Gauss pulses j33[ for example) 
were deduced. 



In Ref. [1^ Besieris et al introduced a novel integral representation for synthesizing 
those LW's. This bidirectional plane wave decomposition is based on a decomposition of 
the solutions of the scalar wave equations into the forward and backward traveling plane 
wave solutions and it has been shown to be a very natural basis for description of LW's (see 
Ref. |3l[ for example). 



The FWM's have been interpreted as being related in a special way to the fie 
moving on a complex trajectory parallel to the real axis of propagation 



d of a source. 



lam^. This 
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observation linked the FWM's with the works by Deschamps |15Cl| and Felsen [151 1 where 
the Gaussian beams have been described as being paraxially equivalent to spherical waves 
with centers at stationary complex locations. 

There has been a considerable effort in finding the LW solutions in other branches of 



physics, spanning various differential equations like spinor wave equation 



hyperbolic systems like cold-plasma equation [23], Klein-Gordon equation 



li 



2211 fist-order 



15|. 



In 1988 Durnin |57] published his paper on so called Bessel beams (see for example 



Ref. |98!l for an earlier publication on the topic). The ide a at t racte d much interest and 

and Bessel-X waves 



the Besse l bea ms and their pulsed counterparts - X- waves Il04l 

r~inn 

|121L 11221 1123| - became the research field of its own rights. In this conte xt the issue of 



109, 



115, 



the superluminal p ropagation of a class of LW's has been considered in Refs 

HQ HQ. 

It has been shown that the FWM's can be described as monochromatic Gaussian beams 
observed in a moving relativistic inertial reference frame 

Propagation of optical pulses or beams without any appreciable drop in the intensity and 
spread over long distances would be highly desirable in many applications. The obvious uses 
could be in fields like optical communication, monitoring, imaging, and femtosecond laser 
spectroscopy, also in laser acceleration of charged particles. Due to this general interest the 
experimental generation FWM's and LW's has been discussed in numerous publications (see 
Refs. H Q 

and references therein). 

The most widely discussed approach has been to use directly the principle of Huygens 

I I 

and launch the LW's from planar sources ^43]. However, it appears that each point source 
in such array must (i) have ultra-wide bandwidth and (ii) be independently drivable as the 
temporal evolution of the LW's generally is of the non-separable nature. Due to the present 
state of the experiment this approach has not been realized even in radio-frequency domain 



(it has been realized in acoustics |4(1 l4l||). 



In an another approach it has been shown that the LW's can be launched by the so-called 
Gaussian dynamic apertures, that are characterized by an effective radius that shrinks from 
an infinite extensions at t ^ — cxd to a finite value at t ^ 0, then expands once more to 



an infinite dimension as t 
Gaussian apertures 



oo 



46 



47 



45[ or by the spectrally depleted (finite excitation time) 



It has been shown, that the field from an infinite line source contains a FWM component 
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17| and the LW's can be generated by a disk source moving "more slowly than the speed 
of light" silsil- 

In optics none of those methods is feasible. A practical general idea for optical gener- 



ation of LW's was described in [53|, |5^ where it was shown that the angular dispersion of 
various Bessel beam generators can be used to produce t he ne cessary coupling between the 
monochromatic components of the LW's. In Refs. j52, 5^, 123 1 we also presented the experi- 



mental evidence of the feasibility of this approach. In particular, in Ref. |j56| we constructed 
an optical setup for generation of two-dimensional FWM's and obtained results from in- 
terferometric measurements of the generated wave field that exhibit all the characteristic 
properties of the FWM's. 

In this review we make an effort to give all the essential results of the field an unified 
description in the way that we present them using the Fourier decomposition methods ex- 
clusively. In doing so we unify the notation and transform the mathematical representation 
where necessary. 

The review is organized as follows: 

In the preliminary Chapter |n] we introduce the necessary integral representations for the 
solutions of the Maxwell's equations and scalar homogeneous wave equation. Predominantly 
we will use the Fourier representation of the free-space wave fields. 

In Chapter IIIII we deduce what in our opinion is the physically most comprehensive 
representation of the FWM's and LW's - we will show, that the necessary and sufficient 
condition for a free-space wave field to be propagation-invariant is that its support of angular 
spectrum of plane waves is of a specific form. Several additional conclusions on the properties 
of the LW's will be drawn. 

In Chapter IIVI we give an outline of the properties of the known (published) LW's. The 
material in this section is important, because, to our best knowledge, this is the first time 
where the optical feasibility of certain well-known closed-form LW's is estimated - we will 
see that majority of the known LW's, including the original FWM's, are not realizable in 
optical domain. 

In Chapter 13 we generalize the theory of the propagation-invariant propagation into the 
domain of partially coherent wave fields - we define the conditions for the propagation- 
invariance of the mutual coherence function of the wideband, stochastic, stationary fields. 
The theory also gives a means of estimating the effect of spatial and temporal coherence of 
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the source light on the properties of generated fields and is used in the analysis of the results 
of our experiments. 

In Chapter IVll we present the general idea of the optical generation of LW's. First 
of all, the setup for the generation of simplest special case - optical Bessel-X pulses - is 
introduced. Then we show that in Fourier picture the optical generation of FWM's can 
be resolved to applying specific angular dispersion to the Bessel-X pulses and discuss on 
the finite energy approximations of the FWM's. Also, the optical generation of partially 
coherent propagation-invariant wave fields is discussed. 

In Chapter I VI II we present the results of the experiments on optical LW's carried on so 
far. In particular, we report on experimental measurements of the whole three-dimensional 
distribution of the field of optical X waves - Bessel-X pulses - and provide the experimental 
verification of the optical feasibility of FWM's. 

In Chapter IVIIII we give an outline of our work on self-imaging pulsed wave fields - it 
appears, that certain discrete superpositions of the FWM's can be used to compose spa- 
tiotemporally self-imaging wave fields that carry non-trivial three-dimensional images. 

II. INTEGRAL REPRESENTATIONS OF FREE-SPACE ELECTROMAGNETIC 
WAVE FIELDS 

In this preliminary chapter we introduce the necessary integral representations for the 
solutions of the homogeneous Maxwell's equations and scalar homogeneous wave equation. 
Only the free-space wave fields are considered, i.e., the wave fields under investigation do 
not have any sources (except perhaps at infinity) and they do not interact with any material 
objects. As we will see, such an approach is suitable for our purposes. 
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A. Solutions of the Maxwell equations in free space 



In SI units the source-free Maxwell equations can be written as 

OH 

V X E = -f^o-g^ (la) 

V X H = (lb) 

V ■ E = (Ic) 

V ■ H = , (Id) 

E and H being the electric and magnetic field vectors respectively, fio is the magnetic 
permittivity of free space, Eq is the electric permittivity of free space. As it is well know, 
in this special case the components of the electric and magnetic field vectors satisfy the 
homogeneous wave equation 

V=-^^)H(r,i) = 0. (2b) 

In Eqs. ()2a|) and ()2b|) only two of the six field variables are independent and the Maxwell 
equations have to used to solve for the other, dependent field components. 

The general solution of the scalar wave equations (j2b| can be expressed as the Fourier 
decomposition as 

E (r, t) = J / duj dk ^ (k, lu) exp [ikr - itut] (3a) 

(27r) J-oD J J J-oo 

H (r, t) = J / dto dkn{k, u) exp [ikr - icot] , (3b) 

(27r) J-oo J J J-oo 

where £ = {Sx,£y,£z) and H = {Hx,'Hy,Hz) are plane wave spectrums of the electric and 
magnetic field. Specifying, for example £y as two solutions of the scalar wave equation 
we get from V ■ E = that 

£z (k, u}) = [k^Sx (k, Lo) + kySy (k, u)] (4) 
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OH 

and from V x E = —Un^^ 

^ dt 

(k, tu) = ^ [KkyS, (k, iu) + - kl) £y (k, u)] (5a) 

ny{\i,u) = — ^ [ ( A;2 - fc2) (k, uo) + fc^fcj^^j, (k, u)] (5b) 
7^2 (k, a;) = -^ky£^ (k, cj) - k^Ey (k, a;) . (5c) 

If we substitute the Eqs. (0} - (j3cj) in ()3a|) and ()3b|) we have a general solution of free-space 
Maxwell equations as a superposition of monochromatic plane waves. 

The other approach is to determine the vector potential A as the solution of the homo- 
geneous wave equation - if we use the Coulo mb g auge and no sources are present the scalar 



potential is zero and the fields are given by 



164| 



E(r,t) = -^A(r,t) (6a) 
B(r,t) = V X A(r,t) (6b) 

Alternatively, we can deterrnine the Hertz vectors 11 from the homogeneous wave equation, 
then the fields are given by |l65l | 

E (r, t) = V (V ■ n(^)) - /.oV X |n(™) - i^n(^) (7a) 
H (r, t) = V (V ■ n(-)) - .oV X |n(^) - ^^n(-). (7b) 

The choice of the vector components of the Hertz vectors and vector potential generally 
determine the polarization properties of the resulting vector field. 



B. Plane wave expansions of scalar wave fields 

If we assume, that the general solution \1/ (r, t) of the scalar homogeneous wave equation 
can be decomposed into the Fourier superposition of plane waves as 



dt I I I dr ^' (r, t) exp [-ikv + iojt] , (9) 

-OO J J J — OO 
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the inverse transform yields 

^'(r,t) = j/ duj dk ^ (k, cu) exp [ikr - iwt] . (10) 

{2tc) J-oo J J J-oo 

The Eq. fllO|) together with the condition 

ki+ki+ki=e=[^y (11) 

which assures, that the Fourier representation satisfies the wave equation (jHJ, is the general 
source-free solution of the scalar homogeneous wave equation that will be used in this review. 
The representation (fTUI) leads to Whittaker an d W e yl typ e pl ane w ave expansions (for the 
discussions on this topic see for example Refs. jlS.S^ [l8fij | and |l^). 



1. Whittaker type plane wave expansion 

The dispersion relation (|TT|l can be inserted into (fTnjl as a delta function 5(A;^ — kl + 
ky + k^) so that the integration over u yields 

or 



where 



" (2^ III "^^"^"^^y"^^' 

X A {kx, ky, kz) exp [i [k^x + kyU + k^z — kct)] . (13) 



A {kx, ky, kz) — 2A;^ ikx^ ky, kz) (1^) 



If we also introduce the cylindrical coordinate system (p, z, (p) in real space and spherical 
coordinate system {k,9,(j)) in /c-space the Eq. (fT^ yields 

\E''(r, t)= / dkk'^ / d0sin6' / A (fc sin cos(/), sin 6* cos</), fccos^) 
Jo Jo Jo 

X exp [ik (x sin 6^ cos (p + y sin 9 sin (p + z cos 6' — ct)] (15) 
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(here and hereafter we omit the normahzing constants in front of the integrals of this type). 
We can also expand the radial dependence of the angular spectra as the Fourier series 



A [k sin 9 cos k sin 6 cos 0, k cos 9) = An {k, 6) exp [ 

n=— oo 

and get another form of p5|) 

^ /"OO i*iv 

\l/ (p, z, if, t) = exp [ibmyj] / d/cA;^ / d6' sin 9 
n=o -^0 -^0 

X An {k, 9) Jn (/cpsin 9) exp [i/c cos 9 — ct)] , 



(16) 



(17) 



where J„ () is the n-th order Bessel function of the first kind and we introduced the polar 
coordinates in real space {p,z,ip), so that \E' (p, z, yj, t) = {p cos (p, p simp, z,t). In the 
radially symmetric case only the term n = is taken into account in Eq. (|17|) and we have 



(18) 



* (p, z, t) = / dkk^ d9sm9 Ao{k,9) 
Jo Jo 

X Jo {kp sin 9) exp [ik {z cos 9 — ct)] . 



If we define x = k sin 9 and again use the Fourier series expansion of the radial dependence 
of the angular spectrum, the representation (fT^ yields 



°° POD POO / 

* (p, z,ip,t)=J2 / dA;, / dxxAn I 

n=oJ~oo Jo V 



X^, arcsin — — 



X 



(19) 



Jn {xp) 6xp [±m<y9] exp i ykzZ — ct^/)^~+~k^ 
Again, in the radially symmetric case only the term n = is taken into account and we have 

/oo POO / 

dkz / dxx^o ( a/^ITx^, arcsin ■ 
-oo Jo \ 



X Jo ixp) exp ilkzZ- ct^ x? + kl 



(20) 



2. Weyl type plane wave expansion 

If we use the dispersion relation (jllj) to eliminate the variable k^ instead, then Eq. (jl(J|) 
can be given the following form 



°° roo roo 

^ (p, z, ip,t) = ^ exp [±in(p] / dk dxx 
n=o -^0 -^0 

X {k, x) Jn ixp) exp \i (z^k^ - x^ - kct 



(21) 
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which is the Weyl type superposition over the plane waves (see for example Ref. |l63l | for a 
thorough treatment). 

The Weyl type spectrum of plane waves is often derived as the Fourier transform of the 
wave field in plane z = 0. In contrary, the Whittaker type superposition is calculated as 
its three-dimensional Fourier transform over the space. Note however, that the distinction 
between the two is not clear for wideband wave fields, as the calculation of Weyl represen- 
tation requires the knowledge of the evolution of the wave field on the z = plane for all 
times [see Eq. Q]. 

C. Bidirectional plane wave decomposition 

The bidirectional plane wave decomposition was introduced by Besieris et al in Ref. 
and it has been proved to be useful for description of LW's. It is based on a decomposition 
of the solutions of the scalar wave equations into the forward and backward traveling plane 
wave solutions, in this representation the general solution to the scalar wave equation can 
be written in the form (Eq. 2.22 of Ref. |l4 |) 



where r] = z + ct and C, = z — ct. Even though the Eq. (j^^ differs noticeably from the 
Fourier decomposition, there is one to one correspondence between these two through the 
change of variables 




X Jn ixp) exp [±mv9] exp [— iaC] 6xp iPr] 5 [ a[3 



(22) 



kz = a — j3 



(23a) 




(23b) 



or inversely 




(24b) 



(24a) 
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Consequently we can write 

oo 

■ip (k, u) oc exp [imyj] 



n=0 



X Cri 



1 /UJ 



Note that the delta function constraint 



4a/3 = 



in Eq. (j22|l in the Fourier picture reduces to 



-) -*^-x' = o. 

C / 



For circularly symmetric wave fields the bidirectional expansions yields 

1 



oo POO 



^ (P, C, ^, V') = -— 3 / da/ d/3 / dxxCo Jo(xp) 

(iTTj Jo Jo Jo 



X exp [— mC] exp ijSi] 



(25) 



(26) 



(27) 



(28) 



III. A PRACTICAL APPROACH TO SCALAR 
FWM'S 



A. Propagation invariance of scalar wave fields 



1. The angular spectrum of plane waves of the FWM's 



First of all, in literature the term FWM has been used mostly with the following closed- 
form solution of the scalar homogeneous wave equation: 



^/ (P, (f, t) = exp [i(3 {z + ct)] 



ai 



Am (ai + i() 



exp 



(29) 



(Eq. (2.1) of Ref. [iSl)- The Weyl and Whittaker type plane wave spectrums of this wave 



field have been derived in Refs. 
reads 

.(/) 



eyi ana 



Ay ^ {k, 0) = - exp 



and, omitting the normalizing constants, the latter 
aik (cos 9 + 11 



6{k~ kcosO -2/3) . 



(30) 
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FIG. 1: On the geometrical interpretation of the parameters /3 and of the supports of angular 
spectrum of plane waves of FWM's (gray line). 



In this respect one can say that the following derivation of the angular spectrum of plane 
waves of the FWM's is nothing but the different interpretation of the results already pub- 
lished. However, the alternate emphasis in the theory, described in this section (and pub- 
lished in Ref. j^), have proved to make the difference if the optical generation of the FWM's 
is under discussion. Also, the term FWM will be redefined in what follows. 

Consider the general solution of the free-space wave equation represented as the Whittaker 
type plane wave decomposition Eq. (fT7j) 



The integral representation of fundamental FWM's can be derived from the condition that 
the superposition of Bessel beams in Eq. pip should form a nondispersing pulse propagating 
along the z axis. In terms of group velocity dispersion of wave packets this condition means 
that the on- axis group velocity = du/dk^ should be constant over the whole spectral 
range. This restriction allows non-trivial solutions only if we assume that the cone angle 
in relation kz = kcos6 is a function of the wave number, i.e., one can write 6{k). The 
corresponding support of the angular spectrum of the plane wave constituents of the pulse, 
i.e., the volume of the fc-space where the angular spectrum of plane waves of the wave field 
is not zero, is a cylindrically symmetric surface in the fc-space and the angular spectrum can 
be expressed by means of Dirac delta function as A„ {k, 6) = (k) 6 [9 — 6 (k)]. 




oo 



n=0 



X An {k, 0) Jn {kp sin 0) exp [ikz cos — iurt] . 



(31) 
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The condition 



= c 



dk 



~~Tr (kcosOp (k)) 
c dk 



dkz 

where constant 7 determines the group velocity, yields 



-1 



(32) 



= 7* - 2/37, (33) 

where the integration constant 2(3 is defined as the wave number of the plane wave component 
propagating perpendicularly to the z axis, i.e.. Op (2/3) = 90° (see Fig. [T]for the geometri- 
cal interpretation of the parameter (3 in fc-space, the choice is consistent with Q| for 
example). Thus, we can write 

cosgp(A:) = ^^^7^^^ (34) 

K 

or 

kp (9) = (35) 

7 — cos 6 

It appears in section UlI A 31 that for a subclass of special cases the above definitions are 
not appropriate as the corresponding supports of the angular spectrum of plane waves do 
not intersect with the k^ axis. Then one should determine an alternate integration constant 
^ from the condition 9f = 180°, this choice yields 

...W J^*-«<^'+^'-''«^° (36) 
[7*-«(7-l), it{<0 

(see Fig. ^for the geometrical interpretation of the parameter ^ in /c-space). Thus, we can 
write 

cos Op [k) = (37) 

rv 

or 

kp (9) = (38) 

7 — cos 

so that 

213 = (39) 

(as we always have /5 > 0). 

The definitions (jH^ or give the angular spectrum of plane waves in Eq. (p?T|l the 
form 

Alf^{k,e)=B,,{k)5[9-9p{k)] (40) 
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and 

) {k, 6) = Bn [Of m 5[k-kF {6)] , (41) 

correspondingly (see Fig. |21of the section ITlI A 31 for the set of special cases). 

As kz = kcos6, our result 1)401) is consistent with the support of angular spectrum of 
plane waves of the original FWM's in Eq. (fHUj). the constant 7 just generalizes to include also 
FWM's of different group velocities. Thus, we can conclude that the physically transparent 
condition indeed determines the support of the angular spectrum of plane waves of the 
FWM's. 



2. Integral expressions for the field of the scalar FWM's 

With the angular spectrum of plane waves ()40|) we can eliminate variable 6 in integral 
dSH) and get 



(p, z, if, t) = exp [±my)] / (ikk'^smOpik) 



n=0 



X Bn (k) Jn [kp sin Op {k)] exp [ik {z cos Op (k) — ct)] , 
using (jH^ we can write 

(p, z, (p, t) = exp [—i2'yPz] exp [±m(y?] / d/c /c^ sin (fc) 



n=0 



X Bn (k) Jn 



kpJl 



k 



exp [ik ijz — ct)]. 



Alternatively, we can eliminate k by means of Eqs. ()^T|) and get 

00 

^> F (p, z, ip, t) = exp [±m(y9] / d^sin^fc^(6') 

n=0 

X Bn [9f (k)] Jn [kp {0) p sin 6] exp [ikp {9) {z cos 9 - ct)] 



again pi)) transform the equation to 



(p, z, (p, t) = exp [±myj] / d^^ 



n=0 



2/37 



X Bn [9f (k)] Jn 



2/97P sin 6* 
7 — cos 9 



exp 



7 — cos 9 
7 — cos 9 



sin 9 

z cos 9 — ct) 



(42) 



(43) 



(44) 



(45) 
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[note that analogous expressions can be written using Eqs. (jHUjl - 

The apphed condition (jH^ imphes that the longitudinal shape of the central peak of the 
pulsed wave field in Eqs. (ji^ - do not spread as it propagates in z axis direction. 
From the integral expressions it is also obvious that the pulse do not spread in transversal 
direction. However, the wave field has what has been called the "local variations" - the 
term exp [— 'j27/9z] in ()43p implies that only the instantaneous intensity of the wave field 
is independent of the propagation distance, in what follows we refer to such wave fields as 
propagation-invariant. 

It is important to note that in Eqs. (jlOl), PT|l and - the frequency spectrum is 
arbitrary. Thus, the necessary and sufficient condition for the propagation-invariance of the 
general pulsed wave field (j31|) is that its support of angular spectrum of plane waves should 
be defined by Eq. (jH^ or (jHTj) . The statement can also be inverted and one can say that the 
wave field is strictly propagation-invariant only if its support of angular spectrum of plane 
waves is defined by Eq. (jl0|) or PT| - indeed, in Eq. p2|l any other choice would lead to 
the group velocity dispersion and the pulse would inevitably spread as it propagates. This 
also implies, that all the possible solutions of scalar homogeneous wave equation that have 
extended depth of propagation as compared to ordinary Gaussian pulses (see next chapter) 
should be considered as certain approximations to the FWM's. 

Now, the closed-form expressions like (j^^ are very convenient in numerical analysis, 
however, limiting ourselves to the set of available closed-form integrals of (j^2|) - P5|) is not 
reasonable by any means. In this review we use the term "focus wave modes" (FWM) for 
all the wave fields that can be represented by the integral expressions (jl^ - whereas 
the closed-form expression ()29|) will be called the original FWM. 

3. A physical classification of FWM's 

The recognition, that the spatiotemporal behavior of FWM's is determined only by the 
support of their angular spectrum of plane waves enables one to give a straightforward 
general classification to the FWM's. 

Note, the dispersion relation 




(46) 
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FIG. 2: The physical classification of the FWM's in terms of sections of the cone + ^2 ~ = 
in (x, kz, k) space. The first two columns depict the sections of the cone from two viewpoints, the 
corresponding supports of angular spectrums of plane waves are depicted in third column. 

can be interpreted as a definition of a cone in (x, kz-, k) space (see Fig. EI). In this context 
the specific supports of the angular spectrum of plane waves of FWM's in Eqs. (jS^ - (|HHj) 
have a geometrical interpretation as being the cone sections of ()46|) along the planes 

kz = -fk- 2/?7 (47) 
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dnSl) or 

K{k)=^k-^{^±1) (48) 

(IHUj) . It can be seen that the possible supports of the angular spectrums of plane waves 
can be divided into four explicit special cases (see Fig. |2I) that can be taken as the natural 
classification of the FWM's: 

1. /? = 0(^ = 0),7<1, the support is a cone in fc-space, typical examples are Bessel-X 
pulse and X-pulse (the case 7 = 1 corresponds to plane wave pulse); 

2. /? 7^ (^ 7^ 0), 7 = 1, the support is a paraboloid in A;-space, typical example is 
FWM's, propagating at velocity of light; 

3. /? 7^ (^ 7^ 0), 7 > 1, the support is an ellipsoid in /c-space, the group velocity of the 
FWM's satisfies < c; 

4. /? 7^ (^ 7^ 0), 7 < 1, the support is hyperboloid in /c-space, the group velocity of the 
FWM's satisfies > c; 

Thus, there is barely four general types of strictly propagation-invariant solutions of the 
scalar wave equation. This point has to be stressed as the straightforward basic idea we set 
forward here is often elusive in the general literature and numerous closed form LW's have 
been set forward. 

a. Pulsed localized wave fields in dispersive media It should be noted at this point that, 
in principle, the approach can be used to derive propagation-invariant wave fields for linear 
dispersive media. In this case we should replace kz by kznioj) in Eqs. (p?^ - (jSSl), n{uj) 
being the refractive index of the medium. This modification yields the following equation 
for the support of the angular spectrum of plane waves of the FWM in linear dispersive 
media: 

V c / n [ijj) ^ 

c being the velocity of light in vacuum. Eq. defines the support of angular spectrum 
of plane waves to the wave field that propagates without any longitudinal or transversal 
spread in linear dispersive media. This approach - to use predetermined angular dispersion 



(49) 



to suppress the longitudina' 
been already used in Refs. 



(arid transver sal) dispersion, though differently formulated, has 



122, 



123 . 



187| for example. 
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4- The temporal evolution of the FWM's in the radial direction 

The temporal evolution of the FWM's in the radial direction can be given a conve- 
nient mathematical interpretation. Namely, as Chavez-Cerda et al noted in Ref. the 
monochromatic Bessel beam can be represented as a superposition of so-called Hankel waves 



(p, z, t) = H^^^ (xp) exp [ik,z - iujt + irrnp] (50a) 

(p, z, t) = H^^ ixp) exp [ik,z - iujt + irrnp] , (50b) 

where 

ixp) = Jrn (XP) + ^Nn. (xp) (51a) 

H"^^ ixp) = Jm ixp) - iNrn ixp) (51b) 

are the m-th order Hankel functions and denotes the m-th order Neumann function (the 
Bessel function of the second kind). For monochromatic wave fields the two solutions define 
the diverging and converging wave in xy plane, in other terms, they form the "sink" and 
"source" pair. In those terms the m-th order Bessel beam can be written as 

Jm ixp) 6xp [ikzZ — iut + irrnp] = (52) 
[H^^ ixp) + {XP)] exp [ik,z - iut + imp] 

- this is a standing wave that arise in the superposition of the two Hankel waves (note how 
the singularity of the Neumann functions at the origin is eliminated). 

This approach can be easily generalized for the wideband wave fields - in this case the 
superposition of the monochromatic Hankel beams form a converging or expanding circular 
pulse in the xy plane. If we also use condition (jiUI) we get the pulse that corresponds to 
the radial evolution of the FWM's. The results of a numerical simulation of its behavior are 
depicted in Fig. andEb- 

Note also, that the radial wave that propagates away from the z axis is generally not 
propagation invariant. Indeed, if we follow the arguments of the section Qll A II for radial 
propagation we can write the condition of propagation-invariance as 
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FIG. 3: (a) Typical spatiotcmporal field distribution of a FWM; (b) The temporal evolution of 
the FWM in radial direction as the superposition of the pulsed Hankel beams (blue solid lines), 
the amplitude of the corresponding carrier-frequency monochromatic Hankel beam is added for 
comparison (red dotted line); (c) The support of the angular spectrum of plane waves of a wave 
field that is propagation- invariant in radial direction (see text). 

where constant 7^ again determines the group velocity. Specifying the integration constant 
^ again from the condition 9p (^) = 180° we can write for the support of angular spectrum 
of plane waves 




(54) 
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Thus, we can write 

sin#(A:) = ^^^^i^^^ (55) 

(note that in this context ^ > 0). 

A typical support of the angular spectrum of plane waves defined by Eq. (jK^ is depicted 
in Fig. inb- So, the FWM is propagation-invariant in both the z axis direction and radial 
direction only in the special case ^ = where we can write 7p = ^1 — 72. This consequence 
will be given a further interpretation in section IIII C 11 



5. The spatial localization of FWM's 

For most practical cases there is no closed-form integrals to Eq. (021) • Consequently, we 
have to deal with integral transforms and the straightforward numerical simulation of any 
realistic situation may be a tedious task (this is especially true for general LW's where the 
double integrals have to be computed). However, for LW's there is a simple method for 
qualitative estimate of the resulting wave fields, based on three-dimensional Fourier trans- 
forms (the monochromatic case of the a pproach Av as introduced by McCutchen in Ref. jissl ] 



and has been used for example in Refs. 

Let us start with Whittaker type plane wave decomposition in Eq. ()13|) and set t = 0: 

I* I* 1*00 

(x, y, 2;, 0) = j j J dKdkydk, 

X A {kx, kykz) exp [i {k^x + kyU + kzz)] . (56) 

Obviously we can write for the field on z axis the relation 

dk, exp [ik^z] I / / dk^kyA (fc^, kyk,) \ , (57) 
00 (.J J —00 ) 



so that 



where 



/oo 
dK exp [ik,z\ A^y {K) , (58) 
-00 

/• /*oo 

A^y{K)= \ / dk.AkyA{k^,kyk,) (59) 
J J —00 



and from the definition of one- dimensional Fourier transform, we can write 



^' (0, 0, z, 0) = 27r [A.,y (fc,)] . (60) 
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FIG. 4: The Fourier transform estimation of the spatial shape of the FWS's (see text). 



Here T^^ [...] denotes the inverse Fourier' transform in fc^-direction and the integral can 
be thought of as the projection of the angular spectrum plane waves onto the z axis (see 
Fig. 13]). Similarly we can write for the filed in xy plane at z = 



^' (x, y,0,0) = / dk^dky exp [ik^x + ikyy] 

J J — oo 

X I /" dk-A{k^.kyk-:}\ , (61) 



so that 



v]/' (x, y, 0, 0) = {27if J=-} [A, (A;., ky)] , (62) 
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where 

/oo 
dkzA {kx, kykz) (63) 
-oo 

and J-'~y [...] denotes the two-dimensional inverse Fourier transform. 

Now, having in mind the table of basic one- and two-dimensional Fourier transforms and 
the general properties of Fourier transforms, the knowledge of the defined projections of 
angular spectrum of plane waves onto the kz axis and k^ky plane allows one immediately 
estimate the general shape of the wave field on z axis and xy plane respectively. If we also 
note that in studies of the propagation-invariant wave fields the estimates are valid over the 
entire z axis (for space-time points 'jz — ct), the approach can prove to be very useful. 
Let us specify the frequency spectrum of the light source s (k) as the Gaussian one: 

1 



s (k) = exp 



(64) 



where ko denote the mean wave number of the wave field and 0"^ is determined from the 
pulse length Tg of the corresponding plane wave pulse as 

cr. 



(65) 



2\/21n2 

From the known character of the angular spectrum of plane waves of the FWM's we can 
approximate for the Gaussian profiles of the kz and x projections of the angular spectrum 
of plane waves 



cosOp (ko) 
(^k 



(66a) 
(66b) 



sm9F{ko) 
respectively. 

The spectral profile of the fc^-projection of the angular spectrum of plane waves then 
reads 



Axy {k^) oc exp 



{kz - kzof 



with the FWHM (full width at half-maximum) 



Akz ~ A/c cos 9f (ko) — 



2V2 In 2 



(Jz 



where 



2 V2 In 2 



(Jk 



(67) 



(68) 



(69) 
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The corresponding intensity profile is 



[A^y (/c^)] oc exp 



2a? 



with FWHM 



(7^2 V2 In 2. 



(70) 



(71) 



For the field in transversal direction we can give a good estimate by recognizing that the 
intensity profile on xy plane has the Bessel profile that is multiplied by an envelope. The 
profile of the latter can be estimated by the ID Fourier transform of the projection of the 
angular spectrum along an axis and we can write 

1 



A-z {kx, ky) oc exp 



-o^p(x-Xo)' 



with FWHM 



Ax ~ AA;sin6'ir (A;o) 
The corresponding intensity profile reads 



2V2 In 2 



T [Az {kx, ky)] oc Jo (kpsmOF (ko)) exp 



(72) 



(73) 



(74) 



with FWHM 



Ap 



Op2\j2 In 2 - the envelope 

- the central peak of the Bessel function 



2x2.405 



(75) 



(see Fig. m for an illustration of the description). Note, that as we can write the ratio 

— = — — = tan6'F (ko) 



Ap 



(76) 



for the pulse widths in the two directions, we at once can deduce that for optically feasible 
FWM's [6p (ko) <^ 1] the central peak is better localized along the z axis. 



B. Few remarks on properties of FWM's 

1. Causality of FWM's 

In several papers it has been noted, that the original FWM's introduced by Brittingham 
and Ziolkowski in Eq. are not exactly causal as they include backward propagating 
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FIG. 5: On the causal and acausal components of the angular spectrum of plane waves of FWM's. 
The striped region denotes the acausal region of the support of the angular spectrum of plane 
waves. 

plane wave components (see Fig. E)) |3| • This fact is due to the specific frequency spectrum 
(IHUl) that leads to the closed-form FWM's (see the overview in following chapter). In the 

n 

consequent publications (see Ref. [[13]]) Shaarawi et. al. demonstrated, that the parameters 
of the spectrum can be chosen so that the predominant part of the energy of the FWM's is 
in forward propagating plane wave components. 

In the context of our approach this problem has to be considered as ill-posed - as all 
the wave fields that share the support of the angular spectrum of plane waves (pn|) are 
propagation-invariant regardless of their frequency spectrum, we can just choose one without 
the acausal components. 



2. FWM's and evanescent waves 



The second topic that is closely related to the backward propaga ting plane wave compo- 
nents of the original FWM's is the one of evanescent waves (isL l2n|. 

From the practical point of view it may seem peculiar to introduce the evanescent waves, 
the intensity of which decays exponentially, in the context of the propagation-invariant wave 
fields where the depth of the propagation usually extend over several meters. However, the 
evanescent waves appear indeed in a Weyl picture of the FWM's. Indeed, from Eqs. fl21|) 
and pHjl one can write 



(77) 
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FIG. 6: The integration contour in Weyl picture of angular spectrum of plane waves. The vertical 
part of the contour where the imaginaty part of the angle 9 is nonzero cancels out in integration. 



for the angular spectrum of plane waves of the FWM's so that the field can be written as 

(78) 



^ /"OO /"OO 

^ (p, z, ip,t) = exp [±inif] / dk dxx B„ (k) 

n=0 



X Jf 



ixp) S ^k — a/A;^ — — 2/? exp i ^za/ k"^ — — ckt^ 



In Eq. (ffH|l . for the ranges x < ^ the integration is over homogeneous plane waves. For 
X > k, the wave vector of the plane waves is purely imaginary and the integration is over 
the evanescent waves Q, 1^^. The situation may be more apparent if we transform to 
variables x^ k ^ 6, k and write (for cylindrically symmetric component only for brevity) 

¥^^{p,z,t)= I dkk"^ I d^cos^sin^ (79) 



or 



Jo Jd± 

X Bq {k) 5[k^ k cos 6 — 2/3] Jq {kp sin 6) exp [±.ikz cos — iuit] 



2/3 



f+w . f 2/?sin^ So , 
(p,z,t) = / d^- ^ 
Jd± 



7 (1 — cos( 



(80) 



X Jn 



2/5psin^ 
7 (1 — cos 6 



exp 



±i- 



2/3 



7 (1 — cos I 



[z COS =F ct) 



Here " +" stands for forward propagating plane wave components and " — " stands for back- 
ward propagating plane wave components and the integration is carried out along the con- 
tours of complex 6 plane, x/k = sin {Or + iOi) (see Fig.lH)). Also, if the analysis is carried 
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out for wave fields the angular spectrum of plane waves of which has forward and backward 
propagating components, the total wave field can be written as 

where subscript denotes homogeneous component of or i.e., < Oji < 2tt, 
9j = and subscript "ev" denotes evanescent components, i.e., 6^ = 7i/2, 9j < 0. It has 

n 

been shown |.18], that for the evanescent components of a free field one has 

= , (82) 

SO that the Weyl forward and backward propagating components add up resulting in the 
source-free solution in Eq. (j42p . 

Again, in our approach the frequency spectrum is chosen so that the wave fields do not 
have any backward propagating components. Consequently, the integration is only along the 
real part of the D^. Also, it is quite clear that for the free-space wave fields the presence of 
the evanescent waves in the integration (jHUj) is rather a peculiarity of the Weyl type angular 
spectrum of plane waves. For example, if we write the Weyl picture of a plane wave pulse 
propagating perpendicularly to z axis, the corresponding Weyl picture obviously do contain 
evanescent components. However, there is no physical content in those components. 



3. Energy content of scalar FWM's 

As already noted, the total energy content of FWM's is infinite 
energy content is calculated as 



/OO poo p27r 

dz / dpp / d(f\'^F{z,p,ip,t)f . 
OO Jo Jo 



Indeed, as the 



(83) 



In the Fourier picture, the Parseval relation and the angular spectrum of plane waves in 
Eq. (PT|) can be used to yield 

2 



POO PIT 

Utot = J2 dk de (k) 5[e-eF (k)] 

POO pTV 2 

= J2 dk de Bn (k) 5^ [e - 9f {k)] 
„ Jo Jo 



so that 



Utot = OO 



(84) 



(85) 
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due to the 5"^ in the integrand (here and hereafter the tilde on angular spectrum indicates that 
the factor fc^ sin 9 is included into the spectrum). Obviously the relation (jH3j) is valid whenever 
there is a delta function in the definition of the angular spectrum of plane waves. Also, it 
has been proved that any wave field that is strictly propagation-invariant has necessarily 
infinite total energy P, Q . 

The second important energetic parameter of the LW's is their energy flow over a cross- 
section per unit time - obviously, any physically feasible wave field has to have a finite energy 
flow. In terms of the previous section and using the two-dimensional Parseval relation this 
quantity can be calculated as 

I* poo 

^.y= / / dKdky\A,{k^,ky)\\ (86) 

J J — oo 

where {k^, ky) is again the projection of the angular spectrum of plane waves onto k^ky 
plane. Obviously the quantity is necessarily infinite, if only the projection of the angular 
spectrum can be written in terms of delta function in k^^ky plane. Otherwise the energy flow 
is finite, provided the function [k^, ky) is square integrable. The comparison of Figs. |21 
and 13 shows that the FWM's generally have finite total energy flow. 

In literature the flnite energy LW's have been constructed for exampl e by means of 
superpositions of FWM's j^lsil and by applying flnite time windows 
section. I VI CI we will describe our approach to this problem as described in Ref. j54| . 

C. Alternate derivations of scalar FWM's 

1. FWM's as cylindrically symmetric superpositions of tilted pulses 

As to demonstrate the efficiency of the integral transform representations in describing 
the properties of FWM's, we give yet another description of FWM's (Ref. js^). 

Let us represent FWM's as the cylindrically symmetric superpositions of the interfering 



pairs of certain tilted pulses (see also Ref 
can be expressed as [see Eqs. (fT^ and (pn|) 



■• 3)- 



In this representation the field of the FWM's 



'^F{p,z,t)= d(t)[T {x,y,z,t;(l)) +T {x,y,z,t;4) + 7i)] (87) 

^0 

d(f)F' {x,y,z,t; 0) , 







33 



where T (x, y, z, t; 0) denotes the field of the tilted plane wave pulses, that in the spectral 
representation are given by 



where A {k, Op (k) , 0) is the angular spectrum of plane waves of the wave field and the 
angular function 9f {k) is defined by Eq. From Eqs. (jHTj) and (jHHj) we get 



An example of a tilted pulse with Gaussian frequency spectrum corresponding to approxi- 
mately Tg ~ 4/s in Eq. (|HH|) is depicted in Fig. [7^, the corresponding superposition of two 
tilted pulses in Eq. dHHl) and FWM in Eq. are depicted in Fig. Eb). 

In this representation the properties of FWM's can be given the following interpretation: 

1. The localized central peak of FWM's is simply the well-known consequence of taking 
the axially symmetric superposition of a harmonic function. Indeed, the interference 
of the two transform-limited tilted pulses in Eq. (j89p gives rise to the harmonic inter- 
ference pattern, the transversal width of which is proportional to the temporal length 
of the tilted pulses (jHHj) . The central peak arises due to the constructive interference 
of the tilted pulses along the optical axis, formally, the cos () function in Eq. (jH^ is 
replaced by Jq () in Eq. ^ [see Fig.[7|D]; 

2. The nondispersing propagation of the optical FWM's wave fields can be given an 
alternate wave-optical interpretation. Namely, it can be seen from Fig. |H1 that in large 
scale the longitudinal length of the tilted pulses depends on the distance from the 
optical axis so that the tilted pulses have a "waist" (this claim is identical to that 
given in section IIII A4I that the radial wave propagating toward the z axis and back 
is not propagation- invariant). The relation (j34|) essentially guarantees, that the waist 
propagates along the optical axis and do not spread - in this case the central peak of 
the corresponding cylindrically symmetric superpositions, FWM's fl42|) . also remains 
transform-limited; 




(88) 



X exp [ik (x cos sin Op {k) -|- y sin sin Op {k) + z cos Op {k) — ct)] , 




(89) 



X cos [k sin 6p (k) {x cos (j) + ysin 0)] exp [ik {z cos 9p {k) — ct)] . 
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FIG. 7: (a) On the field distribution of tilted pulses; (b) comparison of field of the superposition 
of a pair of tilted pulses (in left) and of the corresponding FWM (in right); (c) on the difference of 
phase and group velocities of the FWM's (see text). 

3. The local variations of the central peak of the wave field, noted for example in Ref. 0, 
can be explained as the result of the difference between the phase and group velocities 
along the optical axis - as can be seen from Fig. ^ the pulse and phase fronts of the 
tilted pulses are not parallel; 
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FIG. 8: The large-scale behaviour of the spatial shape of the modulus of the tilted pulses. 



4. The group velocity of the wave field can be set by changing the parameter 7 in Eq. (jH^)- 
The Fig. gives this effect a wave optical interpretation - it can be seen, that the on- 
axis group velocity of the wave field directly depends on the angle between the phase 
front and pulse front and on the direction of the wave vector of the mean frequency. 

It is easy to see, that all the presented arguments are equally valid for the superpositions 
of tilted pulses in Eq. ()89|1 and for its cylindrically symmetric counterparts - FWM's. Thus, 
we can state that the defined interfering pair of tilted pulses possess all the characteristic 
properties of FWM's. In fact, the physics behind the two wave fields is similar to the degree, 
that we will call the wave field ()89|) 



as two-dimensional FWM (2D FWM) in what follows. 

We end this section by noting that the special case of this approach can be used to discuss 
the properties of X-type pulses (Ref. In this case 9f {k) = const = Oq and we have the 

interference of two plane wave pulses: 



00 



F(x,z,t) 



dA; Bq (k) cos [kx sin Op (k)] exp [ik {z cos Op (k) — ct)] 



(90) 




(91) 



X exp [ik {x cos (f) sin Oqi/ sin (p sin 6q + z cos 6q — ct)] 



so that 




00 



dk Bo (k) cos [kx sin 60] exp [ik {z cos 60 — ct)] . 



(92) 



(see Fig. ^ 
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FIG. 9: (a) On the phase and group velocity of a plane wave pulse propagating at angle relative 
to z-axis; (b) comparison of the field of the superposition of a pair of plane wave pvilscs (in left) 
and of their corresponding cylindrically symmetric superposition - Bessel-X pulse (in right); (c) 
on the group velocity of Bessel-X pulses. 
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2. FWM^s as the moving, modulated Gaussian beams 



In literature the closed- 
with the use of the anzatz 



orm expression ()29|) for the original FWM's have been derived 



^ (x, y, z, t) = exp [ipfi] F' {x, y, () , 



(93) 



where fi = z + ct and ( = z — ct. With the wave equation (jS)) reduces to the Schrodinger 
equation for F' 

{A^ + 4i(3dc)F'{x,y,C) = (94) 
which, assuming axial symmetry, has a solution of the form 



F' (P, C) 



exp 



ai + i( 



Ani (ai + iQ 

so that one can write the solution similar to the FWM's in Eq. ()29|1 



^/ (P> P^ = exp 



ai 



exp 



To give the FWM a more convenient form one can use the transform 



1 



1 



— i- 



with which the Eq. 



can be shown to yield 
Wo 



X exp 



4na, (C) 

«I (0 ^(0 V \ai 



where 



ai (C) = Wo 
R{0 = C 



1 + 



ai 



1 + 



and 



If on e co mpares the Eqs. 



(95) 



(96) 



(97) 



(98) 



(99a) 
(99b) 

(100) 



Ref. 



()100|) to those of the monochromatic Gaussian beam (see 
166l | for example) one can see that, the FWM's can be interpreted as moving, modulated 
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FIG. 10: A numerical example of the field of the original FWM's. 



Gaussian beams for which ai {Q and R {Q are the beam width and radius of curvature 
respectively and Wq is the beam waist at = (see Refs. for relevant descriptions). 

Now, several interesting consequences can be drawn at this point. Most importantly, this 
formal analogy between the FWM's and Gaussian beams is very conditional and even mis- 
leading in some respects. First of all, the constant (3 is by no means the carrier wave number 
of the FWM's as one might expect from the corresponding monochromatic expression - in 
the following Chapter IIVI we will see that the convenient choice of parameter for optically 
feasible FWM's with the carrier wave number fco ~ 1 x 10''^ the parameter is of the order 
of magnitude (3 < 100^. Secondly, the requirement of optical feasibility also implies that 
ai ^ 1 (see Sec. IIVB|) and with this condition the original FWM's (see Fig. [TUI) typically 
do not resemble that of the Gaussian beam as they appear in the textbook examples. The 
reason for the "abnormal" behaviour is obvious - with the above conditions the direct anal- 
ogy to the monochromatic case, where (3 = 2it/X, yields for the beam waist in Eq. ()1()()|1 
Wq <^ A. So that we have a limiting case where the waist of the Gaussian beam is much 
less than its wavelength - clearly here the different physical nature of the FWM's show up. 

Next we would like to discuss the claim, often encountered in literature, that the original 
FWM's are carrier free wave fields. First of all, in lights of the general physical considerations 
in section IIII A 51 it should be evident that the non-oscillating shape of the central peak in 
Fig. Uni is a direct consequence of the ultra-wide frequency spectrum of the wave field - 
if the pulse length of the corresponding source plane wave pulse is less than the central 
wavelength, the resulting FWM is effectively an half-cycle pulse and in this condition the 
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FIG. 11: On the character of the Gouy phase shift term in the closed form expression of the 
FWM's: The real part of the original FWM (dotted blue line) is depicted for three z-coordinate 
values together with the modulus (solid green line) and real part (solid black line) of an FWM 
with narrower bandwidth (see text). 

concept of carrier wavelength is rather meaningless of course. However, in above sections 
it was shown that the general FWM's are not confined to the one particular frequency 
spectrum. Correspondingly, we can choose a feasible frequency spectrum and the carrier 
free behaviour of the original FWM's should certainly not be mentioned as the defining 
property of the original FWM's, this is just a mathematical peculiarity of a particular 
integral transform table entry. 
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The issue can be given an alternate description if we note that, using the analogy to the 
monochromatic Gaussian beams the term 



G (z, C) = i (^arctan + 2(3^ (101) 



in the expression Eq. (PS|) could be interpreted as the Gouy phase shift [ |l66l |] of the FWM's. 
In previous section lTlI C ll we described the FWM's as the cylindrically symmetric superposi- 
tions of certain tilted pulses. Now, the original FWM's differ from those, depicted in Fig.[7|D 
only by the ultra-wide bandwidths. In Fig. (fTTj) we have depicted the on-axis spatial evo- 
lution of an FWM as described by Eq. (jHEI) and of one of reasonable bandwidth, calculated 
from the Eq. (021) • The comparison of the two waveforms shows that term exp i arctan 
of the phase term G can be interpreted as the remnants of the sinusoidal waveform, lost 
due to the ultra-wide bandwidth and the term exp [i2(3z\ is added as the monothonically 
growing phase factor that is due to the difference between the group and phase velocities of 
FWM's. The latter term is characteristic to the FWM's only - instead of having a single 
focus with accompanying Gouy phase shift or a "frozen" Gouy phase shift as the X-type 
pulses, FWM's have periodically evolving phase shift term. 

The idea of Gouy phase shift, initially introduced in the Fresnel approximation of the 
diffraction theory of monochromatic focused beams, has attracted a r enewed interest recentl 
in the con text of propagation of subcycle Gaussian pulses (see Refs. 



131 



132, 



;erest recently 
133, HQ 



144L Il4fi| ). We believe, that the simple physical interpretation of the term (jlOlll in the 
context of FWM's, as being the result of the difference between phase and group velocities 
of the wave field, might add to the general understanding of the phenomenon. 



3. FWM's as the Lorentz transforms of focused monochromatic beams 

An interesting interpretation to the FWM's can be given in terms of special theory of 
relativity. Namely, in Ref. |3] Belanger demonstrated that Gaussian monochromatic beams 
appear as FWM's (Gaussian packetlike beams) when observed in an inertial system moving 
at relativistic speeds relative to the focused wave. In this short note we would like to give 
an another mathematical representation to this claim. 

Suppose we take a focused monochromatic wave of the form 

/•TT 

m{p^z,t)= I de K (6) exp [iko{z cos e-ct)] (102) 
Jo 
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with angular spectrum of plane waves 

Ao{k,e) = K{e)6{k-ko). (103) 

If we observe the beam from a moving inertial system, the plane wave components of the field 
suffer from the relativistic Doppler shift. As the result, their wave vectors and frequencies 
transform as described by Lorentz transformations. Specifically, the wave number of the 
wave vector and its longitudinal and transversal components in the inertial frame, moving 
at speed V along the z axis, obey equations 

k' = -fiko{l- picosO) (104a) 

k', = -fiko{cos9- Pi) (104b) 

k'^ = kosinO = kr, (104c) 

k'y = kosinO = ky (104d) 



(Eq. (11.29) of Ref. 16J|) where /cq is the wave number of the wave field in rest frame and 



A {V) = - (105a) 

c 



We can use Eq. fll()4a|l to eliminate 9 from Eq. (jl04b|l to get 



and if we define the parameters as 



7(V) = ^ (107a) 

,3(V.W = ^^. (107b) 
we can write for the z component of the wave vector 

k'^ = -7 {V) k' - 2/3 (F, ko) - 7 iV)^ ■ (108) 

Thus, if the velocity of the moving frame is close to the speed of light, the angular spectrum 
of plane waves of the wave field in the moving frame is the one of the FWM that moves in 
negative direction of z axis - for the FWM's we have in Eq. ()33|) 

k^ = ^k- 2/57 (109) 
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FIG. 12: The support of angular spectrum of plane waves of a focused monochromatic beam as 
seen from inertial reference systems moving at different velocities relative to the rest system of 
the monochromatic beam (fco = 1 x 10''^). Due to the relativistiv Doppler shift the direction of 
propagation and the frequency of the monochromatic components of the focused beam transform 
so that the beam is seen as the FWM in the moving reference system. 



and using both ko and V as parameters we can model every possible support of angular 
spectrum of plane waves of FWM's. In Fig. the support of angular spectrum of plane 
waves of the beam as seen from the moving reference system is depicted for various values 
of the speed V and fixed value for /cq- 

Note that an alternate approach to describe the LW's in terms of generalized Lorentz 
transforms can be found in Ref. js^ - in this work it was shown that the superluminal and 
subluminal Lorentz transformations can be used to derive LW solutions to the scalar wave 
equation by boosting known solutions of the wave equation. 



4- FWM's as a construction of generalized functions in the Fourier domain 

In Refs. jisl . Q| Donnelly and Ziolkowski realized, that various separable and non- 
separable solutions to the wave equation can be constructed in spatial and temporal Fourier 
domain by choosing the Fourier transform of the solution of the differential equation so that, 
when multiplied by the transform of the particular differential operator, it gives zero in the 
sense of generalized functions. In the special case of scalar homogeneous wave equation the 
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corresponding relation reads 



ijj 



^(k,a;) =0, 



(110) 



where (k, uS) Q is (3+1) dimensional Fourier transform of the solution of the wave equation 
(jH}. It can be shown that the function of the general type 



V;(k,u;) = S(x,/5)5 



6 



:iii) 



satisfies flllO|) and yields all the known FWM's (in the sense defined in this review). For 
example the choice 

^ ^'"^^ (112) 



(x,/3) = ^exp 



■ 4/3 



leads to the original FWM's. 

One can notice, that if we eliminate the term x^/^P from the delta functions we get the 
condition (jH^ and thus the Eq. ()111|) is yet another transcription of the support of the 
angular spectrum of plane waves, derived in section IIII A 11 



IV. AN OUTLINE OF SCALAR LOCALIZED WAVES STUDIED IN LITERA- 
TURE SO FAR 

A. Introduction 

Over the years a considerable effort has been made to find closed-form localized solutions 
to the homogeneous scalar wave equations. The main aim of this work is to study the 
feasibility of LW's in optical domain. Without debasing the value of those solutions it 
appears, that this approach often leads to the source schemes that are difficult to realize 
even in radio frequency domain. 

Though there has been several publications that provide an unified approach for the 



description of LW's 



15 



31 



J. I35| ■ to our best knowledge, the optical feasibility of those wave 
fields has not been estimated in literature. Moreover, the analysis of the numerical examples 
that have been published in literature show, that authors have often choose the parameters 
of the LW's so that the frequency spectrum is in the radio frequency domain. 

In our opinion in optical domain the best representation for the analysis is the Whittaker 
type plane wave decomposition. First of all, the mental picture of the Fourier lens that 
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produces the two-dimensional Fourier transform of monochromatic wave field between its 
focal planes is often very useful in modeling the optical setups - we precisely know, how and 
in what approximations the elementary components of the Fourier picture, the plane waves 
and Bessel beams, can be generated. Secondly, the approach of the section Ull A 51 allows us 
easily estimate the spatial shape of the wave fields under the discussion. 

In the following overview we define the term "optically feasible" by two rather obvious 
restrictions: 

1. The frequency spectrum of an optically feasible wave field should be in optical domain; 

2. The plane wave spectrum of an optically feasible wave field should not contain plane 
waves propagating at non-paraxial angles relative to optical axis. 

The latter requirement can be justified by a very simple geometrical estimate, described 
in Fig. 1221^ if the FWM's has to propagate over distances that exceed the diameter of the 
source more than, say, five times, the maximum angle of the plane wave components in the 
wave field has to be less than 5 degrees. 

Note, that the energy content of most of the wave fields discussed in this outline is 
infinite, thus, they are not physically realizable as such. However, as we will see in chapters 
that follow, in optical implementations the finite energy approximations of the LW's follow 
naturally from the finite aperture of the setups and this approximation do not change the 
general properties of the LW's, so that the two conditions for optical feasibility, posed here, 
are also valid for LW's with finite energy content. 

In our numerical examples we try to optimize the parameters of each LW so that (i) 
the frequency spectrum of the wave fields extends from ~ 0.5 x lO'^rad/m (~ llOOnm of 
wavelength) to ~ 2 x lO'^rad/m (~ SOOnm of wavelength), so that ak ~ 3.8 x lO^^m for 
which the length of the corresponding plane wave pulse is ~ 3/s - the shortest possible pulse 
length available, (ii) the plane wave with central wavelength propagate approximately at the 
angle ~ 0.2° relative to the propagation axis, giving o^jop ^ 0.003 for the approximate ratio 
of pulse widths in xy and z direction. Note, that specifying the frequency range and cone 
angle of the Bessel beam of central wavelength completely determines the parameter f3 - for 
7 = 1 we have (3 ~ AOrad/m (clarify section UlI A Ij) . 
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FIG. 13: A numerical example of a FWM with the parameters 7 = 1, /? = 40^, ai = 1.4 x 10~~m: 

(a) The angular spectrum of plane waves in two perspectives; (b) The frequency spectrum of the 
FWM (black line) , the frequency spectrum of an optically feasible wave field (green line) , the angle 
Opik) as the function of the wave number (dashed blue line); (c) The field distribution of the 
FWM. 
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B. The original FWM's 



With Eqs. ()24a|) and ()24b|) the angular spectrum of plane waves of the original scalar 
FWM's in Eq. ^ 



^f{p,z,t)= exp [if3n] 



exp 



Am (ai + iC) ^[01 + 
can be derived from its bidirectional plane wave representation [Q| 



(113) 



TT 



Co(a,/3,x) = 2M/5-2/?)exp [-aa^] 



givmg 



TT 



^0 (k, &) = - exp 



aik (1 + cos( 



5{k-kcose -2(5) 



;il4) 



;il5) 



(see Fig. UnK)- Inserting the angular spectrum (jll5j) into integral representation of the type 
(fT7|l yields 

aik (1 + cos6'i? {k)) 



dk /cexp 



X Jo [A;p sin Op ik)] exp [iA; {z cos 6^^? (/c) — ct)] 



(116) 



(as compared to (fTTj) here we have taken into account the 1/k term that appears in (fT^ as 
to be consistent with for example) so that the frequency spectrum of the superposition 
can be written as 

aik (1 + cos 6'i? (A;)) 



B (k) = fcexp 



;ii7) 



[the significance of the factor ksinOp {k) will be discussed in following sections]. 

The frequency spectrum B (k) in Eq. ()117|1 has two free parameters, oi and f3, the latter 
having the same definition as in Eq. (jHUjl of section IIII A II As we already noted in the 
introduction of this chapter, the choice 7 = 1 together with the frequency range and cone 
angle of the Bessel beam of central wavelength determines j3 = 40rad/m. The single free 
parameter is ai and a single parameter does not allow to approximate for any realistic light 
sources - Fig. [THb shows a typical spectrum that can be modeled in terms of Eq. ()117|1 as 
compared to the optically feasible frequency spectrum specified in the introduction of this 
overview and it can be seen, that the bandwidth of the wave field is far beyond the reach 
of any realistic light source. In fact, due to this extraordinary large bandwidth the original 
FWM's in Eq. ()113|) are essentially half-cycle pulses, as already noted in section IIII C 21 
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As deduced in section UlI C 21 fand also in terms of the section llll A 5|) . the parameter ai 
determines the waist of the wave field - in our numerical example ai = 1.4 x 10~^, so that 
the Eq. ^JW\f gives 

Wo = ~ 6 X 10-^m. 

In literature it has been argued, that the FWM's determined by Eq. (I113p are nonphysical 
as the wave field contain acausal components. In the discussion of Ref. jisf it has been shown 
that the acausality can be eliminated by proper choice of parameters ai and (3 - it has been 
shown that if I3ai < 1, the predominant contribution to the spectrum comes from the plane 
waves moving in positive z axis direction. In Fig IT^ it can be seen, that this is indeed 
the case, however the field is still far from convenient for any optical implementation due to 
ultra-wide bandwidth. 

Note, that various closed-form sub- and superluminal FWM's (7 7^ 1) have been derived 
for example in Ref. Q]. 



C. Bessel-Gauss pulses 



The Bessel-Gauss pulses were introduced by Overfelt in Ref. ||33|| (see also Refs. 
5 lip. In this publication it was shown, that the scalar wave field 







50 



^BG {P,f^X) 



ai 



Kaip 



ai + i( veil + iC 



X exp 



exp 



exp [iPfi] 



A(3 [ai + iC] 



(118) 



ai + i(_ 

where the physical meaning of the parameters ai, (3, also ( and fi is consistent with the 
previous discussion. The expression has an additional free parameter k, as compared to the 
FWM's in ()113|1 . in fact, the latter is the special case of the former in the limiting case 
K — > 0. The Bessel-Gauss pulses were further investigated in Ref. Q| where it was shown, 
that in Fourier picture as in Eq. Q the spatiotemporal Fourier transform of the field can 
be written as 



tpBG (k, uj) = Ebg ix, P) S 



13- — 

^ 4/3 



5 



where 



(x,/3) = -^/o^^J 



exp 



'1^ 



exp 



aiX 
' 4/3 



(119) 



(120) 
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FIG. 14: A numerical example of a Bessel-Gauss pulse optimized for optical generation with the 
parameters a = 40000^, ai = 5 x 10~^m, (3 = 40^, 7 = 1: (a) The angular spectrum of plane 

waves in two perspectives; (b) The frequency spectrum of the Bcsscl-Gauss pulse (black line), the 
frequency spectrum of an optically feasible wave field (green line), the angle Op (k) as the function 
of the wave number (dashed blue line); (c) The spatial field distribution of the pulse. 
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(see section ITlI C 41 for the notation). From Eq. ()119|) it can be seen that the support of the 
plane wave spectrum of the Bessel-Gauss pulses is the same as described by Eq. fjMj) [or 
Eq. (|T7j)]. The change of variables in (jl2Up yields for the frequency spectrum 



B{k) 



aiAn^ f naik sin 6 F (k) 



X exp 



2(3 



^{n' + k'sin'Opik)) 



(121) 



where k > 0, ai > and (3 > 0. 

In the original paper the Bessel-Gauss pulses were introduced as the wave fields that 
are "more highly localized than the fundamental Gaussian solutions because of its extra 
spectral degree of freedom" . The additional free parameter is indeed advantageous, however, 
in our opinion not in the sense proposed in this publication - the spatial localization of any 
wideband free-space wave field is directly proportional to its bandwidth and the latter is 
inappropriately large even for the original FWM's (see Ref. 3] for a related discussion). It 
may be the consequence of this general emphasis of the original paper that it is not generally 
recognized that the extra parameter k in Eqs. flllHj] - 1^121^ gives one the necessary degree 
of freedom to fit an arbitrary bandlimited Gaussian-like spectrum - from Eq. 1)1211) it can be 
seen, that the central frequency and bandwidth of the spectra of the pulse are independently 
adjustable by the parameters n and oi respectively. 

Analogously to the discussion in section ITlI G 21 the Bessel-Gauss pulses can be given the 
form that, in some respect, resembles that of the monochromatic Gaussian beam: 



(p, z, = exp l-zpC] 



ai (0 



^0 



— I- 



Pal (C) R (C) 



X exp 
X exp 



alio 'Amo 

H?aiC Pp^ 



(122) 



4/32a?(C) R{C)) 



arctan | — 
ai 



here again 



ai (C) = 

R{0 = C 



1+^ 



1+1^ 




(123a) 

(123b) 
(123c) 
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The general form of the Bessel-Gauss pulses ()122|) is very advantageous in the sense that 
here we can actually write out its carrier wave number - often this quantity is elusive for 
the wideband wave fields. Indeed, around the point C = 0, along the optical axis (p = 0) 
with (jl23a|) we can write for the z axis component of the carrier wave number 

^o. = /5 + ^-- (124) 
This result is actually quite significant, if we once more remind that in literature the FWM's 



42| for example). In lights 



have often been termed as carrier-free wave fields (see Refs. 
of ()124|) we can conclude that the carrier- free behaviour of the FWM's is indeed caused by 
the integral transform table, not by physical arguments. 

In the numerical example in Fig.EJwe have optimized the parameters of the wave field as 
to match the spectral band specified in the introduction of this section. Again, the parameter 
/? is determined by the bandwidth and the cone angle of the central frequency as described 
above. Thus we got: o = 40000^, ai = 5 x 10~^m, (3 = 40^, 7 = 1. The evaluation of 
the Eq. (fT^ yields ko, = 9.80004 x 10^^ and this result is in good correspondence with 
the numerical simulations. 

In conclusion, the Bessel-Gauss pulses are obviously much more appropriate for modeling 
realistic experimental situations. 



D. X-type wave fields 

The X-type localized wave fields are characterized by that for their angular spectrum of 
plane waves (3 = in Eq. (j34j) [or ^ = in Eq. (j37|) ]. This choice implies, that their support 
of angular spectrum of plane waves is a cone in fc-space (see Fig.|21). Consequently, the phase 
and group velocity of X-type pulses are equal (both necessarily superluminal) and the field 
propagates without any local changes along the optical axis. 



1. Bessel beams 



The Bessel beams j57|-j97| are the simplest special case of the propagation-invariant wave 
fields. Being the exact solutions to the Helmholtz equation in cylindrical coordinates their 
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FIG. 15: A numerical example of a monochromatic Bessel beam with the parameters = 0.223 deg, 
ko ^ Ix 10^: (a) The angular spectrum of its plane waves in two perspectives; (b) The delta-shaped 
frequency spectrum; (c) The spatial field distribution of the beam. 

field reads as 

(p, z, t) — Cn Jji {kpsin Oq) exp [m0] exp [ik {z cos Oq — ct)] (125) 

n 

SO that for the zeroth order Bessel beam we have 

(p, z, t) — Jo (kpsinOo) exp [ik (2;cos^o — ct)] . (126) 
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In the Fourier picture, the zeroth-order Bessel beam is the cyhndrically symmetric super- 
position of the monochromatic plane waves propagating at angles relative to z axis, 
correspondingly, their angular spectrum of plane waves in Eq. (fTH|) reads 

Af\k,e)r^ 5{k-ko)Si9-9o). (127) 

The bidirectional representation of the Bessel beam can be found in Ref. jl^ . 

The properties of Bessel beams have been discussed in many publications both in terms of 



IS nave Deen aiscussea m many pu P lications 
and diffraction theory jT^-jo^l and t 



leir properties 



3: 



where 



angular spectrum of plane waves 

are very well understood today. The interest has been triggered in Refs. 
Durnin et al presented them as " nondiffracting" solutions of the homogeneous scalar wave 
equation - they demonstrated experimentally, that the central maximum of the Bessel beams 
propagates much further than the Rayleigh range predicts. 

Note, that though there has been numerous experiments on Bessel beams, they are not 
realizable in experiment in the exact form (|126p - indeed, the analysis of section HUB 31 
immediately shows, that this wave field has both infinite total energy and energy flow over 
its cross-section. We will discuss this point in what follows. 

In this review the Bessel beams appear as the components of the Fourier decomposition 
in Eqs. (fT3jl - p7|) for example. Later in this review we will refer to their most important 
properties in some detail. At this point we just depict its angular spectrum of plane waves 
with the typical field distribution (see Fig. IT^. 

2. X-pulses 



In 104i | Lu et al demonstrated that the choice 

Al^^\k,e) = B{k)6{e -60) (128) 
in representation p8|l with the frequency spectrum 

B{k) = ^exp[-kao] (129) 
yields the propagation-invariant wave field 

(p sin 6*0) + [ao — i {z cos 6*0 — ct)] 
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FIG. 16: A numerical example of a X-pulse with the parameters Oq = 0.223 deg, 7 = 0.99999: 
(a) The angular spectrum of plane waves in two perspectives; (b) The frequency spectrum of the 
X-pulse (black line), the frequency spectrum of an optically feasible wave field (green line), the 
angle as the (constant) function of the wave number (dashed blue line); (c) The spatial field 
distribution of the pulse. 



(see Ref. |l07l | for the description of higher order X-pulses). From the angular spectrum in 
Eq. ()128|) it can be seen, that the support of angular spectrum of plane waves of the X-pulses 
is a cone in fc-space, i.e., all the plane wave components of the wave field propagate at the 
equal angle from the propagation axis. The frequency spectrum of X-pulses in Eq. fll28|) is 
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uniform (see Fig. IT^ ) - the immediate conclusion of the approach of section IIII A 51 that 
the corresponding field should have exponentially decaying behaviour in both z axis and xy 
plane is confirmed in Fig. ITHh . 



X-wave fields have been further investigated in Refs. 



105, 
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topic have been given an overview and general description in Ref. lll|. We mention here 
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108l |. recently the 



the so called bowtie waves that are generally introduced as the derivatives of the X-waves: 

a™^' (r,t) 



mX 



131) 



The derivatives of X-waves have been shown to possess non-symmetric nature and have 
extended localization along a radial direction. In our terms the physical nature of such 
wave fields can be interpreted by applying the derivation operation on the general angular 
spectrum representation of the free-space scalar wave fields in Eq. (fTH|) . We easily get 



Qm^f (r,t) 



dA; 



d9 (sin I 



\-m+l 



2n 



COS 



dx"' (27r) JO JO 

X A (/c sin 6* COS 0, /c sin 6* COS 0, A; COS 6*) 

X exp [ik {x sin 6 cos (p + y sin 6 cos (p + z cos 6 — iut)] 



(132) 



so that the angular spectrum of plane waves of such wave fields is not cylindrically symmetric, 
correspondingly the wave field is a superposition of higher order monochromatic Bessel 
beams as described by Eq. ()19|) for example. 

Due to the exponential shape of the frequency spectrum the X-waves are not appropriate 
for optical implementation. 



3. Bessel-X pulses 



The Bessel-X pulses were introduced by Saari in Ref. 



12l| as the bandlimited version 



of X-pulses. Their angular spectrum of plane waves can be described as 



4^^) {k,e) = B{k)5{e-eo), 



(133) 



where 



B{k) 



CTk [k 



27r V ko 



exp 



al {k - ko 



(134) 
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FIG. 17: A numerical example of a Bessel-X pulse with the parameters 9o = 0.223 deg, 7 = 0.99999: 
(a) The angular spectrum of plane waves in two perspectives; (b) The frequency spectrum of the 
Bessel-X pulse (black line) as compared to the frequency spectrum of an optically feasible wave 
field (green line), the angle 60 as the (constant) function of the wavel number (dashed blue line); 
(c) The spatial field distribution of the pulse. 
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cTfc being defined in ()65p and ko being the carrier wave number, so that for the field one can 
write 



^BX (p, z,t) 



dk B (k) 



X Jo [kp sin 6*0] exp [—ik {z cos 6^0 — ct)] 



(135) 



The integration in Eq. (fTHj) can be carried out to yield |l21 | 



1 



X exp 



(p^sin^^ + d^) 



Jo \Z (d) pko sin 9] exp [ik^d] 



where 



and 



Z{d) = l + 



id 
koal 



d = z cos 6 — ct. 



(136) 



(137) 



(138) 



From the Eqs. p33|) and p34|l it can be seen that, again, the support of angular spectrum 
of plane waves is a cone in fc-space (see Fig.lTTk). However, unlike the X-pulses, the frequency 
spectrum of Bessel-X pulse is Gaussian and it can be optimized to approximate that of our 
initial conditions. Thus, the Bessel-X pulses are optically feasible in the sense defined in 
this chapter. 



E. Two limiting cases of the propagation-invariance 

1. Pulsed wave fields with infinite group velocity 

Consider the special case 7 = of the support of angular spectrum of plane waves 
that reads 

k, {k)=^ = const. (139) 

From the general definition of group velocity in Eq. (jH^ it is obvious, that for this particular 
case we have f ^ = 00. In what follows we give a physical description to the wave fields that 
have such a peculiar property. 
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FIG. 18: A numerical example of a the wave field with infinite group velocity with the parameters 
^ = 6.7 X lO^m, 7 = oo: (a) The angular spectrum of plane waves in two perspectives; (b) The 
frequency spectrum of the pulse (black line) , the frequency spectrum of an optically feasible wave 
field (green line), the angle Op {k) as the function of the wave number (dashed blue line), (c) The 
spatial field distribution of the pulse; (d) Three snapshots of the temporal evolution of the pulse. 

A closed-form solution of the homogeneous scalar wave equation that obeys ()139p can be 
easily found. The angular spectrum of plane waves of the wave field reads 




(140) 
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The substitution in Whittaker type superposition (fTHj) yields 



where 



so that 



X exp 



Bik) 





l]B{k) 



(p, z, t) = exp [i^z] I dkB (k) Jq (^p V'A?^^^) exp [-ikct] . 

J I I 

If we choose B (k) = const = 1 and use the integral transforms [l| 



:i4i) 



(142) 



(143) 



da; Jo ( 6V — ) cos (xy) = { 



-(3/2_fe2-)-7 gin|^g[(y2_J,2)^j iffe<j,<oo 



and 



dxJn ibVx^ 



jo[uwx~-a-]sm{xy) = 
the integral (|143|) can be evaluated explicitly to yield 

exp^i£,z-i£,y/, 



,if 0<j/<6 

(j/2-62)"5 cos[a(j/2_fe2)2l if b<j/<oo 



(144) 



(145) 



^ exp i^z— c'^t^—p'-^ 



if < tc < p 

if p < tC < CX3 



(146) 



Note, that the special case ^ = yields the cylindrically symmetric superposition of plane 
wave pulses propagating perpendicularly to z axis. 

The support of the angular spectrum of plane waves of the wave field (I14(jp is depicted 
in Fig. llSh.. From the estimates of the spatial localization of LW's in section Ull A 51 we can 
expect the wave field to be localized in transversal direction at t = 0, z = and to be 
uniform along the z axis. Indeed, from the Fig. ITHh it can be seen that at this space-time 
point the wave field ()146|) is an approximation to "light filament" along the optical axis. 

The temporal evolution of the wave field is depicted in Fig. ITHH . One can see, that the 
effect of the infinite group velocity is that the light filament is focused only at a single time 
t = and extends from — oo to oo. As to relate to the conventional wave optics, the temporal 
evolution of the light filament is a close relative to that of the plane wave pulse in a plane 
perpendicular to its wave vector. 
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As the wave field includes both non-optical frequencies and non-paraxial angles it is not 
optically feasible as such. However, it can be shown that a finite energy approximation to 
the light filament can in principle be generated by a cylindrical diffraction grating. 



2. Pulsed wave fields with frequency-independent beamwidth 



In Ref. Campbell et al introduced a wideband wave field that is a superposition of 
the Bessel beams the cone angle of which is chosen so that the condition 

k sin 6 {k) = = const (147) 

is satisfied for entire bandwidth. The condition (jl47|) implies, that the transversal component 
of the wave vector of every plane wave component of the wave field is ao, the corresponding 
wave field was called as the pulsed wave fields with frequency independent beamwidth. The 
angular spectrum of plane waves for such choice can be written as 

(fc, e) = B (k) 6{e- arcsin (^^) ) , (148) 

so that the Whittaker superposition in Eq. (fTHj) yields 







(p, z, t) = Jo (ttop) / dkB (k) exp 

where 



tk \zxll- i^Y ~ct 



(149) 



B (k) = aok'^ B (k) . (150) 

The support of the angular spectrum of plane waves of the wave field Eq. (I149|) is depicted 
in Fig. ^1 (in the numerical example the frequency spectrum B {k) is Gaussian with the 
bandwidth corresponding to ~ 6fs pulse). Using the approach of section IlII A 51 one can 
immediately tell the general spatial shape of such wave fields. Indeed, in this case we have 
a simple special case, where the projection of the angular spectrum of plane waves onto the 
kxky-plsnae is delta-ring, correspondingly, the field in transversal direction at z = 0, t = 
should be of the shape of the Bessel function. As for longitudinal shape, its envelope is 
determined by the bandwidth by it Fourier transform, i.e., we should have a slice of a Bessel 
beam. The numerical simulation in Fig. IT^ shows that this estimate is true. Also, one can 
see that the wave field has generally infinite energy flow. 
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FIG. 19: A numerical example of a the pulsed wave field with frequency independent beamwidth 
with the parameters ao = A x 10^^, 7 ~ 0: (a) The angular spectrum of plane waves in two 
perspectives; (b) The frequency spectrum of the pulse (black line), the frequency spectrum of an 
optically feasible wave field (green line), the cone angle of the component Bessel beams as the 
function of the wave number (dashed blue line); (c) The spatial field distribution of the pulse. 

Comparing the support in Eq. ()148|) to that of the propagation-invariant pulsed wave 
field in Eq. (jH^ and (pn|) one can see, that the wave field ()149|) is not propagation-invariant. 
Consequently, the localized part of the wave field spreads as it propagates. We can also 
suggest the best condition for limited propagation-invariance - the comparison of the support 
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in Fig. IT^ to those of FWM's in Fig. |21 implies that for restricted bandwidths the support 
()147p could be optimized to approximate the " horizontal" part of the ellipsoidal supports of 
the subluminal FWM's (7 > 1). 



F. Physically realizable approximations to FWM's 

As it was explained in section ITll A 51 the presence of the delta function in the support of 
the angular spectrum of plane waves of the free-space scalar wave fields necessarily results 
in infinite total energy content of the wave field. Consequently, for all the above reviewed 
wave fields the total energy content is infinite, 

/oo roc p2n 

dz dpp dip F {p, z, (p,t)f = 00. (151) 
-00 Jo Jo 

Here we proceed by reviewing the approaches used in literature to overcome this difficulty. 
In later chapters we introduce the approach that is especially useful for analyzing optical 
experiments. 



1. Electromagnetic directed- energy pulse trains (EDEPT) 



One approach has been to construct various continuous superpositions of FWM's (111311 

nnnn 

over the parameter j3 (see Refs. p,llJ,ll^lM| and references therein), in this case one writes 



LW 



{z,p,t) 



dpAiP) vl>^ iz,p,t;P) 



ai 



where 



s{z,p,t) 



Ani (ai + iQ Jo 
f3p' 



dp A (/?)exp[s {z,p,t)] 



+ i(3{z + ct) 



(152) 



(153) 



ai + iC, 

A {13) is a weighting function and the subscript LW means "localized wave" . As the supports 
of the angular spectrum of the FWM's for different values of parameter (3 generally do not 
overlap and change smoothly in fc-space, the integration indeed eliminates the delta function 
in the expression for angular spectrum of plane waves (see Fig. QUI). It can be shown 3l| that 
Eq. ()152|) yields finite total energy wave field if only the function A (/?) satisfies condition 



1 1"°° 1 
— / d/?|A(/?)|^-<oo 
2ai Jo /5 



(154) 
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FIG. 20: On the effect of integrating over the parameter (3 on the support of the angular spectrum 
of plane waves of the FWM's: (a) The special case of the "mean" support of the angular spectrum 
of plane waves where 7 = 1 {v^ = c), /3 / 0; (b) The special case where 7 > {v^ < c), again, 



(see Eq. 2.8 of Ref. 3l|), i.e., if only /5~^/^A is square integrable. The LW's of the general 
form ()152|) have been called EDEPT solutions of the scalar wave equation. 

From the discussion of previous chapters it is obvious, that the wave field of the general 
form ()152j) are not strictly propagation- invariant. At first glance it may seem surprising 
because (i) FWM solutions with different values of parameter (3 do travel without any 
spread and (ii) all the FWM's overlap in every space-time point as their group velocities 
are equal. However, the effect can be easily understood if we recollect from section UlI C II 
that the phase velocities of the pulses are different leading to the z axis position dependent 
interference and spread of the superposition of the component pulses (see Ref. jj, 01 for 
alternate proofs of this claim). 

n 

a. Modified power spectrum pulse (MPS) The modified power spectrum pulses 
have been introduced by the following bidirectional plane wave spectrum (see Eq. (3.3) of 
Ref. Isil and Eq. (3.32) of Ref. Q) 



p(p0-b 



,9-1 



(m) 



27rr(g) 





exp 



&ai + (p/3-b) a2 , if/5 > ^ 



, if ^ > /3 > 



(155) 
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FIG. 21: A numerical example of a MPS with the parameters ai = 1.4 x 10 ^m, 02 = 4000m, 
q = lO,p = 0.0001, b = 0.002, (7 = 1, /Sq = 40^): (a) The angular spectrum of plane waves in 
two perspectives; (b) The (3 -distribution in bidirectional plane wave spectrum for a = 0; (c) The 
spatial field distribution of the MPS for t = 0. 

Here 02, b, q and p are new parameters and F denotes the gamma function. Using the 
relations (j24a|) and ()24b|) the corresponding Whittaker type plane wave spectrum can be 
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written as (see also Eq. (3.13a) and (3.13b) of Ref. [3l|) 



27rr(g) 



{MPS) 




(X, ^2 



exp 



(fc+fcz)a 



2 ) '^2 



if A;, < f ^ - ^ 

^ 6 4 p 



(156) 







if t ^ P2<^ 

II «2 > 4 



where = y + ^ind the relation (j^^ has been used. The field function of the MPS's 
is described by equation (see Eqs. (3.34) and (1.4) of Ref. jl^ ) 



MPS 



1 



exp 



47r(ai +<) 



a2 + - 



where 



17]. 



(157) 



(158) 



Att (ai + <) 

The comparison of the bidirectional plane wave spectra of MPS p55|l with that of the 
FWM's ()114|) one can see, that the latter is a special case a2 = 0, g = 1 of the former. 
Consequently, the parameter ai in ()155|) has the same interpretation as in case of FWM's 
- it determines the frequency spectra of the wave field. From (jl56p it is also obvious that 
the parameter 02 determines the width of the (3 distribution and parameter b determines 
the central value of (3. As for parameter q, it can be used to optimize the shape of the f3 
distribution. 

A numerical example of the MPS is depicted in Fig. |^ In this example we tried to 
optimize the parameters so as to satisfy the conditions for optical feasibility as stated in the 
introduction of this overview. From the angular spectrum of plane waves in Fig. I?lk one can 
see, that the MPS's generally have the same inconvenience as FWM's - there no freedom to 
choose the frequency spectrum as to optimize for any convenient light source and they are 
generally half-cycle pulses. 

For an interpretation of MPS's as being the field generated by a combi ned point-like 



144, 



146| 



source and a sink placed at a complex-number coordinate see Refs. 

It is not our aim at this point to study the temporal behaviour of the EDEPT solutions, 
thus, the wave field is calculated only for the time t = 0. 
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FIG. 22: A numerical example of a Splash mode with the parameters ai = 1.4 x 10~^m, a2 = 0.4m, 
q = 16, (7 = 1, Po = 40^): (a) The angular spectrum of plane waves in two perspectives; (b) The 
P -distribution in bidirectional plane wave spectrum for d = 0; (c) The spatial field distribution of 
the Splash mode for t = 0. 



2. Splash pulses 



of Ref. 



Splash pulses p appear if one chooses the bidirectional plane wave spectrum as (Eq. (3.13) 



Cq^^^ (a, f3,x) = -^P'^ ^exp —iaai + (3a2 



(159) 
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One can see, that the bidirectional spectrum is similar in the structure as the one of the MPS 
()155p . Here again the term exp [— aai] can be interpreted as the spectra of the "central" 
FWM and the parameters 02 and q determine the distribution function over the parameter 
(3. The integration in the bidirectional plane wave decomposition (j22j) can be carried out to 
yield (Eq. (3.19) of Ref. Q, Eq. (17) of Ref. Q) 

r(g) 



(a2 - ir]) + 



(160) 



Att (ai + 

The wave field has been called as "splash pulse" in Ref. jsl as for its characteristic spatial 
shape. However, in our numerical example we tried once more to find a set of parameters 
suitable for optical generation. It appeared (see Fig. 1^^ . that in this case the angular 
spectrum of plane waves is very similar to that of the MPS's as in Fig. 1211 

G. Several more LW's 

To date, the literature on LW's and on propagation of ultrashort electromagnetic pulses 
is overwhelming and this overview is by no means complete. Our aim was to demonstrate 
the applicability of our approach on most important special cases. 

As already mentioned, in Ref. js^ Besieris et al derived several closed-form superluminal 
and subluminal LW solutions to the scalar wave equation by "boosting" known solutions of 
other Lorentz invariant equations. 

In section IIIIC4[ we already reviewed the approach of solving the homogeneous scalar 
wave equation and Klein-Gordon equation, introduced in Refs. Q| by Donnelly and Zi- 
olkowski. In those works, they also deduced various closed-form separable and non-separable 
solutions to the wave equation. 

In Ref. js^ Overfelt found a continua of localized wave solutions to the scalar homoge- 
neous wave, damped wave, and Klein-Gordon equations by means of a complex similarity 
transform technic. 

The numerous publications that are involved with ultrashort-pulse solutions of the time- 
dependent paraxial wave equation (e.g. isodiffracting pulses) should be mentioned here (see 
Refs. jl32j |- 146 1 and references therein). 
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H. On the transition to the vector theory 



Even though the scalar theory is often used in description of propagation of electromag- 
netic wave fields, generally the solutions to the Maxwell equations have to be used. However, 
the latter approach is generally much more involved. 

In the context of this review, we investigate the free-space wave fields and mostly use 
the angular spectrum representation of the wave fields. In this context the limitations of 
the scalar theory can be easily formulated - the scalar theory is reasonably accurate if only 
the plane wave components of the wav e fiel d propagate at small (paraxial) angles relative 
to optical axis (see Wolf and Mandel jl63 |. for example). In this review we investigate 
the possibilities of optical generation of FWM's (and LW's), correspondingly, the above 
formulated restriction is satisfied in all practical cases and we can restrict ourselves to scalar 
theory. 

(Of course, this is not the case with the original FWM's and LW's published in literature 
- the plane wave components of those wave fields propagate even perpendicularly to the 
direction of propagation and in their exact description the transition to the vector theory is 
obligatory.) 

The vector theory of FWM' s and LW 's has been formulated and used in several publi- 



|31|, 



lOS 



128 . 
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s a 

by Ziolkowski j^| where he used the Hertz vectors of the form 



13l|. The preferred approach has been the use of the 



cations 

Hertz vectors as formulated in Eqs. (fTajl and (f75|) . One can refer to the theory expounded 



n^'^) = z 

n^™) = z ^F, 



(161a) 
(161b) 



where z is the unit vector along the propagation axis and is the (localized) solution of 
the scalar wave equation. With ()161ap and ()161b|) one get T E or T M field with respect to z 
respectively. A more general treatment can be found in Ref. |l3l| . where the Hertz vectors 
are written as the superpositions of the solutions of scalar wave equation as 



PIS 

n(^) = X ^ a^;^m, + y ^ 6(<^)vl/^ + z ^ c^^m,. 



(162a) 
(162b) 
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The computation of the field components using (|7a|) and ()7d|) . although straightforward, 
results in very complex formulas. 

The intuitive analysis of the effect of the transfer to exact vector theory that is more in 
the spirit of this review can be c arried out in terms of the results that have been published 



on vector Bessel beams in Refs. 



125 



126 
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169|. For example, the result in Ref. |l69j | 



reveals, that for TE and TM fields the vector Bessel beams retain their paraxial-Bessel beam 
nature up to cone angles ~ 14° and this result indeed amply justifies the use of the scalar 
theory in this review. 



1. The derivation of vector FWM's by directly applying the Maxwell's equations 



To finish this chapter we nevertheless advance in some extent the second approach men- 



tioned in Sec. Ill A| where we gave the general expression for the plane wave decomposition 
of the solution of the free-space Maxwell equations. 

To find the vector form for the FWM's as described in Eq. ()42|1 we use the Eqs. ()Ha|l - 
(|Hc|) . In correspondence with Eq. we choose 



£x (k, = {k, (p) 5[kz — k cos 6p {k)] 
Sy (k, Lj) = £y {k, (p) 5[kz — k cos Op {k)] 



and this choice yields from Eq. (jHa|) and (j3 
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(166) 



X exp [ik {x sin 9p {k) cos + y sin Op {k) sin + 72; — ct)] , 
where from Eqs. (jl)) - 

£z {k, (f)) = — tan 6p {k) cos (pE^ {k, (p) + sin (pSy {k, (f. 



(167) 
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and 



c/io cosOp (k) 
(l — sin^ 6f {k) cos^ 0) £y (0, k) 



sin^ Op (k) sin coscpSx (k, 0) + 



TYy (/c,0) 



CyUo cosOp (k) 
+ sin^ 6'i? (fc) sin cos (pEy (k, 
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Hz (k, 0) = sin 9p {k) sin (pE^ {k, 0) — cos (pSy {k, > 
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By expanding 



£i {k, 0) = (/c, n) exp [inc, 

n=—oo 



where 



1 



2tt 



Si {k, n) = — d(j)8i {k, 0) exp [-m0] 



the integration over the can be carried out to yield 

1 r°° 

(r, t) = exp [-2Z/37Z] j 

X exp [m0] Lf {k, p, n) exp [z/c (7^ — ct)] 

1 /""^ 
(r, ^) = 7—2 exp [-2iP^z\ / dA; 

X exp [m0] L]^ (fc, p, n) exp [ik (72; — ct)] . 



where 



(fc, p, n) = k'^ sinOp (k) {k, n) Jn {kp sin Op (k)) 



(170) 
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L^y {k,p,n) = k'^ sin Op (k) £y {k,n) Jnikp sin Op (k)) (176) 

and Lf {k, n, p), (/c, p, n) can be express ed as the hnear combinations of Bessel functions 

169! Il77| for relevant discussions). 



of different order, (see Refs. 
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We also note, that in addition to the TM and TE wave fields azimuthally polarized, 
radially polarized and circularly polarized vector FWM's can be derived ^6^. 
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I. Conclusions. 

The main conclusion of this section are: 

1. At this point it should be clear, that all the possible closed-form FWM's can be 
analyzed in a single framework where the support of the angular spectrum of plane 
waves (j34|) - (j4Up is the only definitive property for propagation-invariance. The 
question of whether an integration over the support has or has not a closed-form 
result is the question of mathematical convenience only. 

2. With a proper choice of parameters some of the closed form FWM's (Bessel-Gauss 
pulses, Bessel-X) are well suited for use as the models for simulating the result of 
optical experiments. In contrary, the LW's we reviewed here - the MPS's, splash 
pulses and the original FWM's - are not feasible in this context. Mostly it is because 
of the ultra-wide bandwidth and non-paraxial angular spectrum content of the pulses. 

3. In our opinion, the procedure of modeling finite-thickness supports for finite energy 
approximations of FWM's reviewed in this section lacks a convenient physical interpre- 
tation and to estimate its practical value this topic has to be addressed in the context 
of a particular launching setup instead. 

V. LOCALIZED WAVES IN THE THEORY OF PARTIALLY COHERENT WAVE 
FIELDS 

Every electromagnetic field in nature has some fluctuations associated with it- even the 
purest laser light is not exactly coherent. However, in optical region the fluctuations are too 
rapid for direct measurement. Their existence can be deduced from suitable experiments 
where the correlation between these fluctuations of field variables at two or more space- 
time points are measured. The second-order coherence theory gives a precise measure of 
those correlations for any two space-time points and formulates the dynamic laws which the 
corresponding correlation functions obey. It provides a unified treatment of all well-known 
interference and polarization phenomena of traditional optics. 

Our interest in introduction of the coherence theory is two-fold. First of all, we can 
generalize the concept of propagation-invariance into the more general class of optical fields 
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where the (ideahstic) fully coherent laser light is but a special case. As the consequence we 
not only obtain a more general view of the subject but also a practical approach towards 
the experimental evidence of validity of the theory presented above. 

In what follows we again confine ourselves to free fields only, i.e., within the spatiotempo- 
ral domain the fields under investigation do not contain sources (except perhaps at infinity) 
and they do not interact with any material objects. 



A. Propagation-invariance in domain of partially coherent fields in second order 
coherence theory 

1. General definitions 

Let us start with some general notion on the subject (see Ref. jl6^). The definitive 
characteristics of a (generally) stochastic wave field in second order coherence theory is 
the ensemble cross-correlation function, often named as mutual coherence function for two 
points that can be defined as 

r (n, r2, ti, ts) = (^*(ri, ti)F(r2, h)), . (177) 

In this equation V(r,t) is the (complex) field of a particular realization of a source-free 
scalar field and angle brackets denote the averaging over the ensemble of realizations, in 
essence the function "measure" the correlation that exists between the light vibrations of 
field V (r, t) at the space-time points (ri, ti) and (r2, ^2) (in terms of statistics the two could 
be called processes). The direct calculation of ensemble averages for such two stochastic 
processes generally requires determination of their joint (two-fold) probability densities p2, 
with which one can write 

r(ri,r2,ti,t2) = j Jv;V2P2[V*,ri,ti;V2r2,t2]dV,*dV2, (178) 

where we have denoted the continuous set of field values in the two space-time points as 

Vi = Vir,,ti), 1 = 1,2. (179) 

The joint probability density is generally unknown, however, there are several special 
cases for which the mutual coherence function can be directly calculated. For example, the 
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two processes, Vi and V2 can be statistically independent of each other, in this case the joint 
probability density is completely separable: 



P2 [V*, ri, h; V2r2, ^2] = Pi {V*, ri, h) pi {V2, V2, 12) , (180) 
so that the mutual coherence function is also separable, giving 

r(ri,r2,ti,t2) = j V*p, {V*,r,,t,)dV* J V^2Pi (V^2, r2, ^2) 01^2 

^{V*(r,,h))jV(r2,t2)),, (181) 

where pi is the first order probability density. The second well-known limiting case is 
the complete correlation (or complete mutual coherence) between the two processes. In this 
case the knowledge of one process completely determines the second process, so that the two 
processes as well as the field ^(r, t) are non-stochastic in nature. The first order probability 
density then have to have the form 

p,{V,T,t) = S{V-V4v,t)), (182) 

where Vc (r,t) is the deterministic function of the space-time point. For the second order 
probability density we have 

P2 [V,*, ri, U; V2, r2, t2] = 5 (V,* - (riM)) S (T4 - K (r2,i2)) (183) 

and for the mutual coherence function of the coherent wave field 

r (ri, r2, ti, ^2) = V; {r,M) K (r2,i2) ■ (184) 

The third special case is the stationary stochastic field. The random process is said to be 
stationary if all the probability densities governing the fluctuations of the field are invariant 
under an arbitrary translation of the origin of time, i.e., if 

p,[V{,r^,t]^p^[V*,r,,h + T] (185) 
P2 [V,*, ri, ti; V2, r2, ^2] = P2 [V*, ri, h + T; V2, V2, t2 + T]. (186) 

In this case the mutual coherence function depends only on the difference of time, so that 
it is often written as 

r (ri, V2, t2, h) = r (ri, V2, t2 - h) . (187) 
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Note, that for the statistically independent processes the conditions ()180p and ()185p imply 
that r (ri,r2,ti,t2) = 0. 

The above definitions are of general nature. Let us now constraint ourselves to the 
stochastic superpositions of solutions of scalar wave fields. In our context we could proceed 
as follows. For scalar wave fields the general stochastic process V{r,t) should satisfy the 
homogeneous scalar wave equation and this obviously implies certain restrictions to the form 
of the corresponding mutual coherence function. In the context of present study the most 
convenient way to introduce them is to represent the field as the Whittaker-type plane wave 
expansion: 

V{r,t)= dk del rf0 a(fc, ^, 0) exp [iA; (rn - ct)] , (188) 
Jo Jo Jo 

where a{k, 6, (p) is the stochastic, properly normalized realization of the angular spectrum 
of the wave field and n = [sin 6 cos (p, sin 6 sin 0, cos 6] is the directional unit vector of the 
plane wave. In this representation the mutual coherence function can be expanded to 

r(ri,r2,ti,t2) = (189) 

p POO p P7V p pItv 

X dkidk2 de^de2 d0id02^(A;i,01,^^l,A;2,02,^2) 



J Jo J JO J Jo 

X exp [-iki (rini - cti)] exp [ik2 (r2n2 - 0^2)] , 

where the A is the (Whittaker type) angular correlation function (cross-angular spectrum 
density), defined as 

A {ki, 01, ^1, k2, 02, 92) = {a*iki, 9i, 0i)a(A;2, 62, 02)) (190) 

and the mutual coherence function obeys the coupled wave equations 

1 (9^ 



Vf-^^jr(ri,r2,ti,t2) = 

V^-^^^r(ri,r2,ti,t2) = 0. (191) 

The angular correlation function ^ () is the ensemble average of the product of the com- 
plex amplitudes a () of corresponding plane waves, i.e., it is the measure of correlation 
between two plane wave components of the field. Note, that the stochasticity in this rep- 
resentation lies in the correlation of amplitudes and phases of the strictly monochromatic 



74 



plane wave components in different realizations of the field. In the coherent limit the phases 
and amplitudes depend on the parameters in non-stochastic manner and do not depend on 
the specific realization. 

In practice the ensemble averages of stationary wave fields are often replaced by the 
corresponding time averages - it can be shown that for certain wave fields the two averages 
are equal, the property is referred to as ergodicity. To this point an apparent contradiction 
can be seen in representation ()188|) - the stochastic part of the definition, the angular 
correlation function a (), do not change in time. Of course, this is primarily the consequences 
of the fact that the temporal evolution of a source- free, forward-propagating, scalar wave 
field is fully determined by its values on a plane. However, to understand how the time 
averaging appears in this representation one has to notice that if the field is stationary and 
the correlations die out sufficiently rapidly as ^2 — ^ cxd and |r2 — ri| —>■ oo, i.e., if the 
field is also ergodic, a single realization of the field can be divided up into sections of shorter 
lengths that are uncorrelated and contain all the statistical information about a realization. 
The time average of the field then is the average over those shorter realizations. Note, that 
in this picture the stochastic angular spectrum (and angular correlation function), being the 
Fourier representation of those shorter sections, depend on time and the time average has 
to be interpreted as an average over the angular spectrums of the shorter sections of the 
realization. 

2. Propagation-invariance in second order coherence theory 

The very nature of electromagnetic wave fields implies that their mutual coherence func- 
tion, the correlation of fluctuations of the fleld at two space-time points, should be a propa- 
gating quantity. Hence, we can also deflne the spatial localization of the mutual coherence 
function and discuss its spread, either transversal or longitudinal. And this is the point 
where one can use the ideas of LW's in coherent theory- we should find the conditions when 
the mutual coherence has some localization quality and when the localization is preserved 
during the propagation in free space. 

The mutual coherence function obeys the coupled wave equations ()191|1 . Thus, the lo- 
calization and propagation of mutual coherence function obeys the same laws as the field 
and in complete analogy to the coherent theory in section IIII A II it can be shown that the 
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instantaneous intensity of a particular stochastic wideband field (essentially a trail of pulses) 
propagate without change along the z axis only if the monochromatic components of the 
field are coupled as 

ap^k, 9, 0) = a{k, 6, (f))6 [6 - Op [k)] , (192) 

where the function Op (k) is defined in Eq. El correspondingly, for the longitudinal compo- 
nent of the wave vectors of plane waves we have 

k, = k cos 9f (k) = -fk- 2Pj (193) 

and the field can be expressed as 

/'2n POO 

F'{r, t) = exp [-i2/3-fz] / d0 / dk a{k, Op (k) , 0) (194) 

Jo Jo 

X exp [ik (x cos sin 9p {k) + ?/ sin sin Op {k) + 72 — ct)] , 
As already noted above, the cylindrically symmetric case of the field reads 

/•oo 

F (p, z, t) = exp [-i2jpz] / dk 

Jo 

X b {k) Jo [kp sin 9p {k)] exp [ik (72; — ct)] . (195) 

The Eqs. ()194|) and ()195|) are regarded the definition of FWM's in second order coherence 
theory in what follows. One can see the intimate relevance between the mathematical 
definitions of coherent and partially coherent theory, the only difference being the stochastic 
nature of the angular spectrum of plane waves for the latter. The direct correspondence 
between the two also allows us to transfer the results on finite energy content LW's described 
above, so that the FWM ()194j) or ()195p is essentially the infinite energy content limit of the 
broader class of fields, LW's. Here we proceed by studying the consequences of the stochastic 
nature of the angular spectrum of plane waves of the partially coherent FWM's. 
The angular correlation function for the partially coherent FWM's reads 

{a*p{kueu(pi)apik2,e2,(p2)) = (196) 

Ap {k,, k2, 01, 02) S [6^ - dp (k^)] 6 [02 - Op {k2)] , 

i.e., the cross-angular spectrum vanishes unless 61 = 62 = 6p {k). The general expression for 
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the mutual coherence function of the propagation-invariant partially coherent FWM's reads 

poo poo p27r p27v 

TF{ri,r2,ti,t2) = exp[-i2(3-f{z2 - zi)] / dkidk2 / d0id02 

Jo Jo Jo Jo 

X Af ih, k2, 01, 02) exp {ik2 [r^2n±2 ik2) + (7^2 - ^2)]} 

X exp {-iki [r^in_Li (ki) + (7^1 - cti)]} , (197) 

where r±j = {xj, yj), and n_|_j (kj) = (cos (f>j sin Op (k) , sin 0^ sin 9p (k)), j = 1, 2. The mutual 
coherence function (jl97|) depends on the longitudinal coordinate and time instant through 
expression z'j — ct (or, equivalently, z — v^t) and through the z dependent overall phase and 
can be expressed as 

Tp (ri, r2, ti, ^2) = exp [-22/^7 {z2 - zi)] 

X G (xi, 2/1, Zi7 - cti, X2, ?/2, 2:27 - Ct2) • (198) 

This is the direct analog of the corresponding expression in the coherent theory (jU^ . 

Without the loss of generality we can introduce a partitioning of the angular correlation 
function as 

Af {ki, k2, 01, 02) = V* (A;i, 0i) V (A;2, 02) C (ki, k2, 0i, 02) • (199) 

Since the modified angular correlation function C depends on all variables in non-factored 
manner, this expression still describes the most general case of the non-stationary non- 
homogeneous partially coherent FWM's. The factors V(..) in Eq. ()199|) can be given the 
interpretation as being the square root of the angular spectrum density S {k, 0) multiplied 
by a phase constant, i.e., 

V {k, 0) = {a*{k, Of (k) , (f))a{k, Of (k) , 0))^ exp [iv {k, 0)] 

= [S {k, 0)] ^ exp [iv {k, 0)] , (200) 

where the phase factor v {k, 0) is a real quantity and is essentially the extracted phase of 
the factor V (..). 

As the factors V(..) do not depend on specific realization, the statistical, correlation 
properties of the field have to be determined by the function C , for example, the special 
case where the function is constant corresponds to full correlation- obviously the integrand 
in Eq. ()197|) is nonzero only for the pairs of plane waves that correlate in different realizations 
for the ensemble. 
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To summarize this section we note, that the defining property of the partially coherent 
FWM's is the same as for their fully coherent counterparts: the support of their angular 
spectrum of plane waves obeys the linearity condition in Eq. (|193|) . In following section we 
discuss several limiting special cases of the angular correlation function in Eq. (jl99p . 

Note, that the following approach can be considered as the generaliza t ion o f the discussion 
on monochromatic propagation-invariant wave fields in Refs. 
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B. Special cases of partially coherent FWM's 



1. Coherent limit 

For fully coherent wave fields the phases and amplitudes of the plane waves in Eqs. ()198|l 
and ()199|1 are invariant of the specific realization of the field and all the plane wave compo- 
nents are fully correlated, hence, the non-factored part of the angular correlation function 
can be expressed as 

C {ki, k2, (pi, (P2) const = 1 (201) 
the angular correlation function A factorizes to 

Af (ki, k2, (j)2) = V* (fci, (f)i) V ik2, (202) 

and the mutual coherence function factorizes as 

Tp (n, r2, ti, t2) = F* (n, ti) F (r2, ts) , (203) 

where 

r2Tr roc 

F (r, t) = exp [-i2Pjz] / dk 

Jo Jo 

X V (A;, 0) exp [ik (r^n^ (k) + {-fz - ct))] . (204) 

Thus, the factor V {k, (p) as defined in Eq. p99|) appears as the plane wave spectrum of the 
coherent wave field the multiplier exp [ih {k, 0)] defining the phase of the monochromatic 
components. The field (j2U4j) could be generated from mode-locked femtosecond laser pulses 
in a stable linear-optical setup. 

To conclude this section we note, that the Eq. p97|) implies that the mutual coherence 
function of a nonstationary field is generally time-dependent quantity. However in optical 
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domain the function of the field are far too rapid for direct measurement, so, in optical 
experiments a time-averaged mutual coherence function appears 

(r (n, r2, ti, h)) = ^ dt V* (ri, t - ti) V (rs, t - t^) . (205) 

If we confine ourselves to cylindrically symmetric wave fields the time- averaging can be 
carried out to yield 

/■oo 

Tf (ri±,r2±,A;z, 7 - cr) = exp [-i/37A2] / dA;|V(A;)|^ (206) 

Jo 

X Jo {kpi sin 6p (k)) Jq {kp2 sin Op (k)) exp [ik {jAz — cr)] . 
2. FWM's with frequency noncorrelation 

Consider the special case where the plane waves of different wave number are not corre- 
lated, however, the field is spatially fully coherent at some particular frequency throughout 
a volume, i.e., the field is completely coherent in space-fequency domain. In this case the 
nonfactored part of the angular correlation function Ap has the form 

C(fci,fc2,0i,02)oc5(A;i-A;2), (207) 

so that 

Ap {ki, k2, 01, 02) = V* {ki, 0i) V {k2, 02) S {ki - k2) (208) 

and the field is stationary, i.e., its statistical prop ertie s do not depend on time origin. The 
cross-spectral density of the field also factorized ([l6^; (4.5.73)) 

W (n, r2, k) = W (n, k) U (r2, k) , (209) 

where 

lA (r, fc) = exp / d0V (fc, 0) exp [zfc (rj_n^ [k) + (7^ — ct))] 

00 

= exp [—2/572;] exp [±my9] V„ (/c) Jn \k r^nj_ (fc)] exp [ifc (72; — ct)] (210) 

n=0 

is the temporal Fourier transform of the field and the quantity V (A;) can be found from the 
relation 

/ d0V (/c, 0) exp [iA;rj_n^] = V„ (A;) exp [±my9] J„ (A; rj_n^) . (211) 

n=0 
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So that the cross-spectral density function can be expressed as 

oo 

W{xi, zi,X2, Z2-, k) = exp [-i2P'y {z2 - zi)] ^ exp [±20 (m - n)] (212) 

n=0 

X \Vn{k)f Jn {kpi Sin Of {k)) Jm {kp2 ^vnOp {k)) exp [ik'^ {z2 - zi)]. 

The mutual coherence function can be written either by taking the Fourier' transform of 
the cross-spectral density (j212p 

r,. (ri,r2,r) = (V* (n, t) V (rs, t + r)) 

/■oo 

= / dkW {Yi,k)U{Y2,k)e-'''^\ (213) 
Jo 

or by inserting the angular correlation function into the general Eq. ()197|) . As the result we 
get 

(ri±, r2±, A2, 7A2; - cr) = exp [-i(3'^Az\ 

nliT p2n POO 

X / d0id02 / dA;V*(0i,fc)V(02,fc) (214) 
Jo Jo Jo 

X exp [ik (r_Lin_Li (A^i) - r±2n±2 (^2) + l^z - cr)] . 
The integration over the azimuthal angle gives: 

00 00 

Tf (ri±, r2±, Az, jAz - cr) = exp [-i(3jAz] Yl ^^15) 

n=0 m=0 

poo 

X exp [±i {rupi - mv?2)] / dfcV* (k) Vm {k) 

Jo 

X J„ (kpi sin 9f (k)) Jm {kp2 sin 9f (k)) exp [ik {•yAz — cr)] , 

so that for cylindrically symmetric fields we have 

Tf (ri±, r2±, Az, 7A2 - cr) = exp [-z/?7A2] / d/c | Vo (/c) | 

X Jo (kpi sin 6^^ (fc)) Jq (fcp2 sin (fc)) exp [zfc (7A2; — cr)] , (216) 

this result is identical to the time-averaged mutual coherence function in (j2U6j) . 

In general the full coherence of the light in space-frequency domain can be achieved by 
filtering the incoherent light from a nearly blackbody source by a small (delta) pinhole. The 
corresponding partially coherent FWM's, described by the Eq. p97p are characterized by 
the two properties: 
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1. Their mutual coherence function depends only on the difference Az, i.e., its is invariant 
of the origin on the z axis. This also means, that the mutual coherence function 
propagates without change. 

2. In general the intensity of the stationary field can be expressed by 

J(r) = r(r,r,0) (217) 

then from Eq. ()215|) it can be seen that the averaged intensity of the FWM's is de- 
scribed by 

oo oo 

/ (r) = ^ ^ exp \±i (nv?i - m^<,)\ (218) 

n=0 m=0 

X / dfcV* (/c) Vm (/c) -'n (fcpsin^i;' (/c)) Jm (/;;psin6'F (/c)) 

and this expression do not depend on location on z axis. In particular, if the field is 
cylindrically symmetric we have 

/"OO 

/(r)=/ dk\VQ{k)f .Jl{kp^\^Qp{k)). (219) 

For wideband fields the integration yields a localized on-axis spot, for quasi- mo nochro 
matic wave fields the intensity in near axis volume is a close approximation to that of 
the zeroth-order partially coherent Bessel beam / (r) ~ Jq [kpsinOp (k)). 

3. FWM's with angular noncorrelation 

Consider the special case where the plane waves of different azimuthal angle are not 
correlated (are directionally 6 correlated) but for every particular direction the components 
add up coherently. In this case the nonfactored part of the angular correlation function Af 
of the partially coherent FWM's has the form 

C(fci,A;2,0i,02)oc5(0i-02), (220) 

so that 

Af {ku k2, 01, 02) = V* {ku 0i) V {k2, 02) 6 (01 - 02) (221) 
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and for the mutual coherence function we have 



r^;' (ri, r2, h, t2) = exp [-i2f3'y {z^ - 

nllT POO POO 

X / d0 / / dA;idA;2V*(A;i,0)V(A;2,0) 



^0 ^0 JO 

X exp [-iki [r^in_L (h) - (72:2 - 0^2)]] 

X exp [ik2 [r±2n± (ki) + - ^2)]] , (222) 
The integration over the azimuthal angle gives: 

Tp (ri, r2, ti, 12) = exp [-i2/?7 (22 - ^i)] 

°° POD POD 

X J2 exp [±m0] / / dfcidfcsV: (fci) V„ (^2) 
n=o "^0 

X J„ (/cip2sin6'F (ki)) J„ (fc2P2 sin 6*^7' (^2)) 

X exp [iki (7Z1 - cti)] exp [-ik2 {'yz2 - 0^2)] , (223) 
The expression (I222p can be rewritten as 

p2n 

Tf (n, r2, ti, ta) = / d0^* (0, n, h) ^ (0, r2, ta) , (224) 

where 

POO 

^ {(j), Vj, tj) = exp [-i2p-fZj] / dkV {k, (j)) (225) 

Jo 

X exp [—ik (r_|_jn_|_ (/c) + (7^ — ctj))] . 

In this special case the plane wave constituents for every particular direction add up 
coherently resulting in a pulsed tilted plane wave constituents \1/ (0, r, t) that have a (gen- 
erally) specific profile for every direction 0. The integration in Eq. (j222|) generally yields 
an incoherent angular mixture of fully coherent plane wave pulses and F factorizes only 
asymptotically in the far field. Such fields might be formed from mode-locked femtosecond 
laser pulses in an optical set-up containing elements or parameters which fluctuate (slowly 
as compared to the pulse repetition rate) - moving scattering media, etc. 

The FWM's of this type have the following properties: 

1. Again, the mutual coherence function propagates without any change. However, the 
time origin and the origin on the longitudinal axis is important now as the field is not 
stationary. 
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2. The time averaged intensity can be found from the mutual coherence functions as 

TF{r,r,t,t)= d0 / / dA;idA;2V*(A;i>i)V(A;2,02) X (226) 
Jo Jo Jo 

X exp [-i {k2n± {k2) - kin± {ki)) r± + {k2 - ki) (72; - ct)] , 

so that 

n2'K poo 

7(r) = / d(/) / dk\Vik,(f))f (227) 
Jo Jo 

and due to the noncorrelation of the tilted pulses the recording system see uniform 
intensity distribution. However, the instantaneous intensity V*V in this case strongly 
depends on location on z axis and time origin. 

4. FWM's with angular and frequency noncorrelation 

We conclude this discussion with the special case where all the plane waves are uncorre- 
lated. The non-factored part of the angular correlation function then can be expressed as 

C (A;i, k2, 01, 02) <xS {cl>, - 02) S {k, - k2) , (228) 

so that the angular correlation function Af does not factorize neither with respect of the 
angular variables nor of the frequency: 

Af (ki, k2, 01, 02) = V* (A;i, 0i) V (^2, 02) S (0i - 02) S {k, - ^2) . (229) 

The corresponding mutual coherence function can be expressed as 

(ri, r2, h, t2) = Tf (ri - r2, h - ^2) (230) 

poo r27r 

^exp[-i2(3j{z2-zi)] dk d0|V(A;,0)|' 

JO Jo 

X exp [-ik [n± (k) {r±i - r±i) - {-fAz - cAt)]] , 
The integration over the azimuthal angle gives: 

(ri, r2, h, t2) = exp [-^2/^7 {z2 - zi)] (231) 

00 poo 

x^exp[±m0] / dk\Vn{k)f 

n=0 

X Jn [k {p2 - pi)smeF (A;)]exp{iA; [7(2^2 - zi) - 0(^2 - ^i)]}- 
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The above equations show, that in case of superposition of completely uncorrelated plane 
waves the mutual coherence function F is homogeneous and stationary, i.e., it depends only 
on differences of its arguments. 

The FWM's of this type have the following properties: 

1. Despite of the complete noncorrelation the mutual coherence function of this type 

of partially coherent FWM's can be localized and it propagates without any change. 
Particularly, the field is propagation-invariant in the sense that F depends on the 
longitudinal coordinate z through the difference zi — z\ only. 

2. The time-averaged instantaneous intensity is described by 



i.e., it is uniform in space. 

This is the field, that may be viewed as the opposite of an coherent FWM in variable- 
spatial-coherence optics: its intensity is uniform along all spatial directions and in time, 
whereas the sharply peaked behavior, characteristic to FWM's, reveals itself in the mutual 
coherence functions Vp. 

C. Conclusions 

In this chapter, we have generalized the concept of propagation-invariance into the domain 
of partially coherent wave fields. In the case of partially coherent LW's the propagation- 
invariant, spatially localized field variable is its mutual coherence function. Nevertheless 
the mathematical description is similar to that used for coherent wave fields - the angular 
correlation function of the partially coherent LW's, being the counterpart of the angular 
spectrum of plane waves in the second order coherence theory, is defined using principles 
known from coherent theory. 




(232) 
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VI. OPTICAL GENERATION OF LW'S 



A. Introduction 

Despite the extensive theoretical work carried out on FWM's and LW's, for the long 
time there was very few experimental verification of the concept of propagation-invariance 
of pulsed wave fields. The only experimental verification on LW's was the launching the 
"acoustic directed energy pulse trains" for ultrasonic waves in water 40|,|4l|. Theoretically 
the problem has been addressed in numerous publications [s^] - 5l| verv seriously - in 
a letter to the author of this review of 17.12.2002 Pierre Hillion wrote: "Please do accept 
my apology but it makes so a long time that I am working with focus wave modes that 
I am skeptical on the possibility of man-made focus wave modes; perhaps only Nature in 
cosmic events?" . In the lights of the discussion of the overview in Chapter IIVI one has 
to agree with this opinion, but with the following concretization. In optical domain this 
opinion is true because (i) today we do not have a coherent light source that generates 
half-cycle pulses in visible region and (ii) due to the non-paraxial angular spectrum of plane 
waves the original FWM's are more like the modes of a cylindrical resonator. However, 
in the preceding chapters we have shown, that the two requirements are nothing but the 
peculiarities of a closed-form integral of the corresponding general expressions that describe 
propagation-invariant pulsed wave fields. In this chapter, we demonstrate that the wave 
fields that are essentially the limited-bandwidth modifications of the original FWM's can 
be generated in optical domain. 

Majority of the publications on generation of FWM's and LW's (except of course those 
that are the part of this review) have discussed possibilities of launching them from an array 
(matrix) of (point) sources. This approach has proved to be very involved - as it has been 
noted in literature j,4^, the LW solutions generally cannot be written as the product of a 
function only of time and a function only of space, i.e., the LW solutions are mathematically 
nonseparable in the space-time coordinates. As a consequence, if one tries to use the principle 
of Huygens and launch the LW's from a planar source, it appears that each point of this 
source have to be driven independently. In other words, the nonseparability of the LW's 
implies that the frequency spectrum of every point-source is different and one has to drive 
a separately addressable array of wideband Huygens sources with a function identical to a 
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FIG. 23: The principal scheme of a setup that generates optical FWM's. Here L is a lens and 
AS is a chromatic angular slit. The angular dispersion arises as the consequence of the frequency 
dependence of diameter of the chromatic angular slit. 

LW. The analysis presented in literature demonstrates that if we could build that kind 
of matrix of sources, the launched wave field would be a perfectly causal FWM. The finite 
aperture of the source or the time-limiting of the excitation do not destroy the localized 
propagation of the generated pulses. 

In optical domain this approach has a fundamental drawback - the frequencies of the wave 
fields are of the order of magnitude lO^^Hz and the idea of driving a matrix of independent 
ultra- wideband sources in this frequency range cannot be realized in experiment. The minor 
deficiency is that in the Huygens representation the propagation properties of the LW's are 
not physically transparent in the degree they are in angular spectrum representation. 

As an example of an alternative approach, characteristic in optics, consider a setup con- 
sisting of a Fourier lens of focal length / and of a circularly sjTiimetric mask with transfer 
function 

t{R,k) = s{k)6[R- ftanepik)], (233) 

where R is the radial coordinate on the ring and 6f (k) is the angular function of a FWM 
defined by Eq. ()34j) - this is essentially an annular ring mask, the diameter of which depends 
on the wave number of the light as R{k) = /tan^i? (k) (see Fig. EBI)- It is well known, that 
each monochromatic point-source of wave number k at some radial distance p on the focal 
plane of an ideal the lens results in an apertured plane wave that subtends the angle Op (k) 
relative to z axis behind the lens. The angular spectrum of plane waves of the total field 
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then can be described as (here the infinite aperture is assumed) 



Ao {k, 9) = s (k) exp [iv (k)] 6[9-eF (k)] , (234) 

where v {k) is some phase factor specific to the setup. The support of the angular spectrum 
of plane waves of the wave field is the one of the FWM's pUj) . If we also compensate for 
the phase factor in Eq. ()234p by applying the conjugated phase chirp to the input pulse the 
Fourier representation directly suggest that a transform-limited FWM can be generated by 
illuminating a "chromatic" annular ring mask by a specifically chirped plane wave pulse. 
Though such chromatic mask is not very practical either, the advantage of this approach 



over the Gaussian aperture 49j in optical domain is obvious - this model is essentially static. 



B. Feasible approach to optical generation FWM's 

Let us introduce the general idea in terms of 2D FWM's as defined in section IIII C II 
Consider the pair of plane wave pulses propagating at angles 9o and —6*0 relative to the 
optical axis (see Fig. 124b). Obviously we can introduce a tilt into their angular spectrum of 
plane waves by means of the angularly dispersive elements like diffraction gratings or prisms 
(wedges) and provided that the resulting tilted pulses overlap, we can observe the interference 
of two tilted pulses in near-axis conical volume [see the striped region in Fig. 1^^]. If the 
introduced tilt of the plane wave components is such that the condition in the Eq. ()34|1 . i.e., 

cosg^(fc) = ^^^~^^^ (235) 

rC 

is satisfied for every k and if we ignore for the while the diffractional edge effect that appear 
due to the finite aperture of the optical elements, the interference pattern in the near-axis 
volume is that of the 2D FWM. 

The cylindrically symmetric case (FWM) can be considered similarly - one has to find the 
means to generate the superposition of Bessel beams, the support of the angular spectrum 
of which satisfies condition (j235|) . Also, the use of the angular dispersion of the known 
generators of Bessel beams, i.e., the devices that transform a monochromatic plane wave 
into a Bessel beam, could be an appropriate idea. 

As the first step we should generate the cylindrically symmetric counterpart of the pair of 
plane wave pulses - the Bessel-X pulse (see section IIV D 3|) . The general approach for this is 
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FIG. 24: The general idea for optical generation of the 3D (2D) FWM's: (a) The FWM generator 
consist of an axicon A and of circular diffractional grating G (or two wedges and linear diffractional 
grating respectively if we generate 2D FWM's), the FWM can be observed in the conical volume 
(striped region) behind the diffraction grating; (b) The support of the angular spectrum of plane 
waves of the initial wave field on the FWM generator (solid line). Here and hereafter the dotted 
line denotes the support of angular spectrum of plane waves of the (2D) FWM under discussion, 
the dashed lines denote the bands of the frequency spectrum of the light used in our experiments 
(note the difference in scales between and kz axis); (c) The support of the angular spectrum 
of plane waves behind the axicon; (d) The support of the angular spectrum of plane waves at the 
exit of the FWM generator as compared with the theory. 
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FIG. 25: The Bessel-X pulse generator. Here L stands for lens and AS for angular slit. 



quite straightforward - according to p33|) we have to use the Bessel beam generator without 
any angular dispersion and illuminate it with a wideband pulse. The obvious choice is the 
annular slit in the back focal plane of a Fourier lens (see Fig. 1^^ . The simple geometrical 
considerations show that in this case the wave field behind the lens is the cylindrically 
symmetric superposition of plane wave pulses - for the cone angle of the Bessel beam of the 
wave number k we find sin 6' {k) = sin (arctan {D/ f)) = sinOo, where D is the diameter of 
the angular slit and / is the focal length of the lens, so that the generated wave field can be 
approximated by 



/•oo 

^Bxip,z,t)= / dk s (k) Jq {kp sinOo) exp [ik {z COS 9o — ct)] , (236) 
Jo 

where s (k) is the frequency spectrum of the source field. 

Using the Bessel-X pulses as input, the problem of optical generation of FWM's reduces to 
the modeling of a set of diffractive elements that transform the conical support of the angular 
spectrum of plane waves so that the angular dispersion of the output pulse approximates 
the one described by the condition in Eq. ()235|1 and to the "compression" of the resulting 
wave field by compensating for the relative phases between its monochromatic components 
so that a transform-limited pulse appears on the optical axis. 

There are various Bessel beam generators described in literature such as axicons, circular 
diffraction gratings, etc. (see, e.g., Ref. 57 1 - 6^ and references therein). By illuminating 
those elements with a plane-wave pulse, we obviously get a superposition of Bessel beams, 
the support of the angular spectrum of which is determined by the wavelength dispersion 
of cone angle of the optical element 6{k). For example, an axicon is characterized by 
the complex transmission function exp [ik tan a {1 — n (k))], where a is the angle formed 
by the conical surface with a fiat surface and n {k) is the refractive index of the axicon 
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(see, e.g., Ref. 



74|). The method of stationary phase 



163[, apphed as in Ref. 



0, 



yields 



sin 9 (k) = tan a[l — n (k)]. Hence, the angular spectrum of plane waves of the corresponding 
polychromatic wave field can be approximately described by 



Aq (fc, 9) = s {k) 5[9 — arcsin (tan a (1 — n (/c)))] . 
Substitution into Eq. (fT7|) yields 

/■oo 

^(p, z,t)= / dk s{k) J Q[kp {tan a {I - n{k)))\ 
Jo 

X exp [ik [z cos (arcsin (tan a{l — n (k)))) — ct]] 



(237) 



(238) 



and we can conclude that ideal Bessel beam generators can be used to design wave fields 
with cylindrically symmetric, two-dimensional supports of angular spectrum. Likewise, a 
circular grating yields for the cone angle sin 6* (k) = 27T/kd, where d is the grating constant 
(first-order diffraction is assumed). The polychromatic wave field can be written as 



^'(p,z,t) 



dk s {k) Jo 



'0 

X exp 



ik 



cos ^ 



(239) 



arcsm 



27r 
kd 



ct 



(note, that the Eqs. (j238|) and (j239|) do not count correctly for the longitudinal intensity 
change of Bessel beams behind axicons and circular diffraction gratings [3, 1^ , they are 
brought about just to exemplify the use of the angular dispersion in our problem). 

The support of angular spectrum of a FWM ()34|) cannot be approximated by a single 
diffractive element. However, we will show below that aside from the conventional configu- 
ration, where axicons and circular gratings are used to transform a plane wave into a Bessel 
beam, those elements can also be used to change the cone angle of Bessel beams. This 
property allows us to use a set of those Bessel beam generators as to design more complex 
supports of the angular spectrum. 

The generic features of the diffraction of Bessel beams on circularly symmetric optical 
elements can be understood by means of the followin g mo del. A zeroth-order Bessel function 
in the off-axis region can be well approximated by |l59l | 



Jo (a) 



TT 



+ exp 



TT 

-I \ a 

4 



(240) 
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FIG. 26: A qualitative description of a Bessel beam as a superposition of diverging (dotted lines) 
and converging (solid lines) conical waves in plane z = : (a) free-space evolution of the conical 
components; (b) diffraction on a circular opaque mask; (c) refraction on an axicon. 



Substitution in 



yields 



2 V Tckp sm6 



exp 



7r 



ikpsinO — i — 
. 4 



+ exp 



-ikp sinO + i— j exp [—ikct] 



(241) 



for the Bessel beam in plane z = 0. In this expansion a Bessel beam in this plane is a 
superposition of two conical waves, as shown in Fig. OUk (see also Ref. jfi^]). The evolution 
of those components can be qualitatively analyzed by means of ray-tracing, if we also rec- 
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FIG. 27: (a) A circular diffraction grating on the surface of an axicon - a composite optical element, 
which can be used for generation of FWM's; (b) An optical setup for generation of FWM's. A 
plane wave pulse is incident upon an annular slit (AS). Bessel-X pulse with cone angle 6 behind a 
Fourier lens (L) is incident upon the composite optical element (AG). 



ognize that the characteristic field distribution of a Bessel beam arises in the volume where 
the conical waves interfere. Typical examples, the free-space evolution and diffraction on 
circularly symmetric mask, are depicted in Fig. OUk and Fig. I^bjj . As one can see, such well- 



known features of Bessel beams as finite propagation length j57[ and regeneration behind 
an opaque screen appear immediately. 

Let us place an axicon (or circular grating) onto plane z = (see Fig. HEt). A qualitative 
analysis by means of ray-tracing immediately reveals, that the element changes the cone 
angle of both conical waves. The " converging" conical wave forms a Bessel beam in a conical 
near axis volume, but the cone angle is now different. The "diverging" conical component 
leaves the near axis region. Hence, axicons and diffraction gratings change the cone angle 
of a Bessel beam in the sense that the regenerated Bessel beam after the element is of a 
different cone angle. 

Let us find the cone angle of the resulting Bessel beam for a composite optical element- 
a circular diffraction grating on the surface of an axicon (see Fig. I7fk). We assume that a 
Bessel-X pulse with cone angle ^o, generated by means of an annular slit and Fourier lens, 
is incident upon it (see Fig. ITTb). Snell's law and the grating equation yield the following 
equation for the cone angle of the resulting Bessel beam: 

1 

n {k) 

Here d is the grating constant and n (k) is the refractive index of the axicon material. Sign 



27r 

sin Og {k) = — + n (k) sin ^ —a + arcsin 

rhd 



sin {9o + a) 



(242) 



conventions are chosen so that the angles a, 9o, Oq {k) are positive in Fig. ITTk. First-order 
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diffraction is assumed. If the angles a, are small, so that sinx ~ arcsina; ~ x, Eq. 
yields 

27r 



(243) 



Similar results can be obtained by means of the method of stationary phase. Below we 
use the approach in Ref . [3] • Let the initial field on the diffractive element be the converging 
conical component of a Bessel beam |a/2/ {irkp sin 6'o) exp [ikp sin 60]. The amplitude trans- 
mission function of the element is exp [ikp {2ti /kd + tana (1 — n {k)))]. The corresponding 
Fresnel diffraction integral can be given as 



^ (P, z) = —e 
iXz 



1 ik{z+er 



D 



where 



and 



dp7 (p') exp [ikp {p')] 
'kpp'^ 



J2 



p(p) = Yz-P 



2 V T^kp' sin^o 
27r 



27rJo 



kd 



+ tana (1 — n {k)) + sin^^c 



The first derivative of p {p') yields exactly one critical point p'^, namely 
p'c = z ^— + tan a {1 — n{k)) + sin 6^0^ = z sin 9sp (k) 



(244) 



(245) 



(246) 



(247) 



and, according to Ref. [7^, the leading contribution to (j244|) behaves as 



iXz 



exp UA; ( 2; + 



2z 



D 



(ip'f (P ) exp [ikp (p')] 

/ (p^) exp [^fc/i (p'J] 



(248) 



iXz 



exp 



^k{z + ^^ 



where p^^^ {p'J denotes the value of the second derivative p (p^) at the critical point. Sub- 
stituting ()247j) into ()248|) and omitting some position-independent factors we get 



\E' (p, z) ^ — — exp 

%Xz 



z sin 9sp {k) 



X 



2 y irkz sin (fc) sin 6q 
z sin 6'5p (fc))^ 



27rJn 



/cpj; sin ^^^p (A;) 



X exp 



ik 



2z 



— z sin^ 9sp (k) 



(249) 
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FIG. 28: The experimental demonstration of the use of the angularly dispersive Bessel beam 
generators for change of the cone angle of monochromatic Bessel beams. In the setup: AS, annular 
slit; L, lens; A, "hollow" axicon. In the left and right pane the CCD image of the Bessel beam at 
the plane A and B are depicted respectively. 



so that 

ik— 
2z 



\1/ (p, z) ~ Jo {kp sin 9sp (k)) exp 



exp [ikz cos 9sp (k)] , (250) 



where relation 1 — sin 9/2 ~ y 1 — sin 9 = cos9 is used. We also have k {p /2z) <^ tx 12 for 
the far-field in near-axis region and Eq. (j249|) reads 

^ (p, z, t) ^ Jo {kp sin 9sp (k)) exp [ikz cos 9sp (k)] . (251) 

Thus, the method of stationary phase yields a similar result as the qualitative ray tracing 
analysis- the output field of an axicon (or circular grating), illuminated by a Bessel beam, 
is a Bessel beam. If we assume a small angle limit i.e. sin a; ~ tanx ~ x in Eq. ()247|) . the 
generated cone angle 6'5'p is also exactly the same as predicted by the ray tracing analysis 
in Eq. ()24H|1 and we have shown that the Bessel beam generators can be used to change the 
cone angle of Bessel beams. We also carried out experimental proof of this principle, the 
typical results are depicted in Fig. |2H1 
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The wavelength dispersion of the cone angle of the Bessel beam constituents of a FWM 
Op {k) is determined by Eq. ()35|) . Combining this condition with the dispersion of the cone 
angle of the discussed setup 6^^ {k) (j243p . we get 

arccos = ^ + a (1 - n {k)) + ^o- (252) 

In following chapters we will see that the three free parameters (a, d) can be adjusted so 
that the relation ()252|) is satisfied in a very good approximation in a limited spectral range. 

As to finish this section we note that care must be taken when specifying the frequency 
spectrum B (k) or B (k) = P sinO (k) in Eqs. P?7j) and ^M)- In terms of the 

frequency spectrum of the light source s (k) - generally the choice the transform depends 
on the particular setup under consideration. If we set s{k) = B{k), we integrate over the 
monochromatic Bessel beams with amplitudes s{k). However, for example in the principal 
setup in Fig. EHl it is quite obvious that the amplitudes of the Bessel beams depend on 
the diameter of the annular slit - the larger diameter means larger area of the slit and 
more transmitted energy. Consequently, the "chromatic" annular ring mask has a transfer 
function that is not constant but can be approximated by ~ A;sin6' (k) = x and one should 
use the weighting function B {k) = k sin 6 {k) s (k) in the superposition over the Bessel beams 
or 

Ao {k, e) = k sin 9 (k) s (k) exp [iv (k)] SIO-Of (k)] (253) 

in place of ()234|1 for the angular spectrum of plane waves (here v (k) stands for the phase 
distortions of the setup). As for the special cases of annular slit with constant diameter, 
axicon and circular diffraction grating the analogous arguments show that the replacement 
s{k) = B (k) applies, (note also how the otherwise exponential frequency spectrum of the 
FWM's in Eq. ()115|1 behaves as a Gaussian-like in ()117|1 due to the k sin 6f {k) term in the 
expression of the Whittaker type plane wave expansion ()17j) ). Of course, the generated wave 
field is propagation-invariant regardless of the weighting function in the superposition and 
generally in any optical experiment only the bandwidth of the source have to be concerned 
about. 
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C. Finite energy approximations to FWM's 



We already noted in the overview in Chapter IIVI that the total energy content of FWM's 
is infinite and as such they are not realizable in any physical experiment. We also refer- 
enced some of the most commonly used approaches to derive finite energy approximations 
to FWM's. Obviously, in experimental situation a natural choice is to consider the ap- 
proximations that correspond to the specific launching mechanism. In other words, the 
approximations of the type ()152|) should be given a physical content, in terms of limitations 
of a real experimental setup. In what follows we derive the finite energy approximation 
of FWM's that is due to the finite aperture of the optical system introduced in previous 
chapter (Ref. 3)- 

In Fourier picture the starting point for this discussion is obvious. In this picture any 
realistic (finite-aperture) optical system generates a superposition of apertured monochro- 
matic Bessel beams so that the field in the exit plane the FWM generator can be described 
by (for brevity we restrict ourselves to cylindrically symmetric case) 

(p, 0, 0;k) = t (p) Jo [kp sin Op (k)] , (254) 

where t (p) is the complex-amplitude transmission function of the aperture of the setup. 
One just has to show that (i) the resulting superposition of the apertured Bessel beams 
still represents a wave field that propagate as a FWM, and (ii) the superposition has finite 
energy content. 



1. Apertured (finite energy flow approximations to) Bessel beams 

The simplest mathematical model for the monochromatic Bessel beams is the exact so- 
lution of the scalar homogeneous wave equation in Eq. ()126|1 

B (p, z, t) = Jo {kp sin 6^0) exp [ik [z cos ^o ~ ct)] . (255) 

In this case the beam can be described in Fourier' picture by the single delta-ring in fc-space 
and the transformation of the beam in free space and in optical elements can be easily 
estimated by simple geometrical constructions. However, such closed form Bessel beams 
are not square integrable. The most apparent approach to deduce physically realizable 
approximations to Bessel beams has been to apply a finite aperture to the system. The 
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FIG. 29: On the finite aperture approximations to the Bessel beams. The optical setup is the 
simple Bessel beam generator consisting of annular ring AS of diameter D, and the lens of aperture 
D and focal length /. The striped region denotes the near axis volume where the generated wave 
field approximates closely the behaviour of the finite aperture Bessel beam. 



problem has been discussed in many works, mostly in terms of Fresnel approximation of 
scalar diffraction theory (see Refs. 

and references therein). 

As an alternative, we could use the Gaussian windowing profile and consider the paraxial 
wave equation - this combination yields a closed mathematical form of the beams, that have 
been called the Bessel-Gauss beams (see Refs. and references therein). 

However, the paraxial approximation lacks the physical transparency of the angular spec- 
trum representation. Still, it should be noted, that the Bessel-Gauss beams can be used to 
construct wave fields that behave like LW's propagating in free space. 

We will use the generally accepted facts that (i) applying finite aperture to a Bessel beam 
provides us with a finite energy flow wave field, that is a very good approximation of the 
infinite-aperture Bessel beams fl254j) in a certain finite-depth, near axis volume (see, e.g., 
Ref. 0|) and (ii) that the polychromatic superpositions of those apertured Bessel beams 
approximate very closely the superpositions of " non-apertured" Bessel beams in this volume 
(67l |. Such a behavior can be easily explained in the Fourier picture where a monochromatic 
Bessel beam is a cylindrically symmetric superposition of plane waves, that propagate at 
angle 6 relative to z axis. Indeed, as the apertured plane waves approximate their infinite 
aperture counterparts very closely in their central parts, one can also observe a very good 
approximation to the infinite-aperture Bessel beam in this near-axis volume (see striped 
region in Fig. 1221 and Ref. 8^ for a more detailed description in terms of diffraction theory). 
If the cone angle of a Bessel beam is small, as it is always the case in paraxial optical 
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systems, the apertured Bessel beam would behave as its infinite-aperture counterpart 
for several meters of propagation. 

The situation in Fig. |2nican be modeled by applying the transmission function t (p) of the 
aperture to the Bessel beams. In Weyl picture this operation is equivalent to calculating the 
two-dimensional Fourier transform of the transversal amplitude distribution t (p) Jq{xoP)j 
where xo stands for radial projection of the wave vector of the Bessel beam. Given the Weyl 
type angular spectrum of plane waves of the infinite-aperture Bessel beam 

A{x)=K5ix-Xo), (256) 



where i^' is a constant, the Fourier transform of ()254|) in point z = 0, t = can be found to 
be 

K 

AapB (x) = -—2 T{x)*5{x- Xo) 
(27r) 



^ d0r (^^ x' + xl- 2xXo COS (0 - 0o)) , (257) 



where T (x) is the two-dimensional Fourier transform of the transmission function t (p) and 
* denotes convolution operation (see also Ref. ^^). The argument of the function T () 
in ()257|1 has an interpretation as being the distance between the points (x, 0) and (xc^o)- 
As for all convenient apertures the function T (x) is well-localized around zero, the major 
contribution to the integral ()257|) comes from small values of and one can write in good 
approximation 

Aa,b (X) ^ ix - Xo) . (258) 

The interpretation of the expression ()258|1 is straightforward: the finite aperture gives the 
support of angular spectrum of a monochromatic Bessel beam a finite "width". Exact form 
of the support is determined by the complex-amplitude transmission function, however, the 
well-known set of fundamental Fourier transform pairs gives a good idea of what the support 
of angular spectrum looks like, without any calculations. 

2. Apertured FWM's 

First of all, it has been demonstrated both numerically and experimentally that the super- 
position of apertured Bessel beam behaves like the superposition of their infinite-aperture 
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counterpart in near-axis volume as defined in Fig. 123 Thus we can claim that the sub- 
stitution of infinite-aperture Bessel beams in (pn|) by their apertured counterparts should 
generate a finite energy flow wave field that is a good approximation of the FWM (j42p in 
this finite volume. In other words, in this near-axis volume the apertured wave field can be 
well approximated by the formulas of infinite-aperture wave fields. 

Obviously, the finite aperture has a similar effect on the angular spectrum support of a 
FWM - the delta function in Eq. (jlU]) is substituted by a weighting function and the angular 
spectrum of plane waves of apertured FWM's can be written as 

XF jk) B jk) r'- 

AapF [k, X) = -—2 / 

{2tt) Jo 



X T [ J x^ + Xf {kf - 2xxf (k) cos<^ 



^^^^^^TW-x.(k)l (259) 

where xf (k) = k sin 9p {k). Consequently, the Weyl type plane wave expansion of the wave 
field behind the aperture reads 

'^ApF{p,z,t) = dk dxxAApF{k,x) 
Jo Jo 

X Jo (px) exp ik {z\fW^^^ — ct^ . (260) 
Alternatively, the transformation x = k sin 6 gives the expression ()260|) the following form 

^ApF{P)Z,t)= I dkk"^ / d^^sin6'cos6' A^pi? (/c, A;sin6') 
Jo Jo 

X Jo {kp sin 9) exp [ik {z cos 9 — ct)] . (261) 



The support of the angular spectrum of plane waves ()259j) of the derived wave field is 
depicted on Fig. IHUk . In correspondence with Eq. ()258|) it has a finite "thickness". We can 
outline the main difference between the support of angular spectrum of LW's, proposed in 
this section (I^HT|l and that of EDEPT's (ITH^ discussed in section IIV F 11 The comparison 
of their supports of angular spectrum of plane waves on Fig. IHUk and IHUb (respectively) 
shows, that the transversal "width" in k^ky plane of angular spectrum of plane waves of 
our apertured FWM's is a constant- such a result is a consequence of applying aperture 
with a wavelength- independent complex-amplitude transmission function. On the other 
hand, the transversal width is generally not constant for the superpositions of FWM's in 
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FIG. 30: The comparison of the supports of angular spectrums of plane waves of (a) an apertured 
FWM and (b) of an EDEPT type superposition of FWM's as described in above sections. 



Eq. ()152|) . The description in this section gives the property a straightforward interpretation: 
the corresponding aperture has a wavelength dependent complex-amphtude transmission 
function. In the context of our discussion, where the main goal is optical feasibility, such an 
approach should be regarded as an impractical one. 

The spatial localization of the result ()260p can be estimated by the approach in sec- 
tion nil A 51 - obviously the resulting wave field still has the characteristic narrow central 
peak. Also, in correspondence with the note in the beginning of current section, the field 
of apertured and non-apertured FWM's do not differ noticeably in the near-axis volume at 
t = 0. Thus, if the aperture of the system is reasonably large (several millimeters) the only 
qualitative effect of the finite aperture is the reduced propagation length of the wave field. 

There is several ways to estimate the propagation length of the apertured FWM's. The 
simplest possibility is still use the geometrical construction in Fig. 123 from which we can 
write 



where fcmin denotes the minimal wave vector of the spectrum of the wave field (note, that 
for the supports of FWM's this generally corresponds to maximum angle max [6p {k)], i.e., 
to the minimal propagation length of the corresponding Bessel beams [clarify Fig. 0^ ]. The 
second possibility is to estimate the propagation length of the apertured monochromatic 
Bessel beams by means of the approach in the section llll A SI fin present case where the wave 
field is es sentially monochromatic, this is actually the McCutchen theory as introduced in 
Ref. llSS^ ). From the Eq. (ffH|) we at once get 
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i.e., we get identical result as in Eq. ()262|) . The third possibility is to estimate the maximum 
group- and phase velocity uncertainty of the finite width support of the angular spectrum 
of plane waves (|257|) . 

In conclusion, in our setup the finite energy (physically realizable) approximations to 
FWM's are introduced quite plainly by the finite aperture of the optical system. In Fourier 
domain this corresponds to smoothing of the delta function in the angular spectrum of 
plane waves of the wave fields. Also, the double integrals in Eq. ()260|) can in some detail be 
handled without excessive numerical calculations. 



3. On finite time window excitation of the FWM's 

In literature, there has been several works on generation of FWM's where the finite total 
ene.,y ha. been aehieved by .,„„t,.g the exeitation t„„e of tKe so.ce a„.ay QH HQ. 



It has been shown, that, indeed, such approach do not destroy the localized propagation 
of the generated pulses. It is intuitively clear, that the finite time excitation results in a 
superposition of the longitudinal "fragments" of the monochromatic Bessel beams. In this 
picture the reasonable time window indeed do not corrupt the behaviour of the central part 
of the wave field, and it is still a good approximation that of the exact monochromatic 
Bessel beam. What the time window do is that it broadens the frequency spectrum of the 
points of the source matrix. Correspondingly, the support of the angular spectrum of the 
generated wave field would look much like the one in Fig. IHUk . Actually, the excitation 
time of our setup is also finite and the two broadening effects - due the aperture and the 
finite excitation time - appear simultaneously, however, the one that is caused by the finite 
aperture is several orders of magnitude stronger. 



D. Optical generation of partially coherent LW's 

The direct comparison of the equations describing the angular spectrum of plane waves 
of coherent and partially coherent LW's in Eqs. - (jlHj) and ()192|1 - ()195j) respectively 
imphes, that the main part of the generating setup should be the same for both cases. 
Indeed, the partially coherent quasi-monochromatic plane waves transform in linear optical 
systems similarly as the coherent monochromatic plane waves, at least the description of the 
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previous chapter concerning refraction and diffraction in grating hold for both cases. The 
quahtative difference hes in the correlations between the plane wave components. 



1. The light source 

According to our general idea of the optical generation of FWM's in Sec. IVI Bl we need 
a well-directed wideband partially coherent plane wave as the initial field in our setup. In 
mathematical limit the angular spectrum of plane waves of such initial field is described by 

a{(t),e,k) = s{k)5[e], (264) 

giving for the field 

1 r°° 

V{r,t) = — dk s{k)exp[ik{z-ct)], (265) 
2vr Jo 

where s (k) is generally stochastic function. Such field is fully coherent in transversal direc- 
tion. In longitudinal direction the coher ence time Tc is determined by the bandwidth of the 
light Ak and the reciprocity inequality 163l | 

TcAk > ^. (266) 

For the ensemble average, we can write 

{s*{kMk2))=S{k)5{k,-k,), (267) 

where S {k) denotes the spectral density (power spectrum) of the light source so that the 
mutual coherence function reads as 

/"OO 

r (Az, r) = / dkS (k) exp [ik {Az - cr)] . (268) 
Jo 

In other words, the mutual coherence function of light field behaves like a plane wave pulse 
of the duration in coherent optics. 

The traditional approach to generate such field is to use a thermal light source, for exam- 
ple the superhigh-pressure Xe-arc lamp, and spatially filter the light by means of a pinhole. 
The well-known drawback of the choice is the hugely reduced signal level as compared to 
laser sources. The alternative is to use a directional white-light continuum source. 
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FIG. 31: The principal setup for generation of FWM's with angular noncorrelation: AS, annular 
slit; LI, L2, L3, lenses; AG, annular grism; D, diffuser; 

2. FWM's with frequency noncorrelation 

As it was explained in section |VB2] (clarify Eqs. (IM^ and ((2101)), the FWM's with 
frequency noncorrelation differ from their fully coherent counterparts only by the lack of 
correlation between the fluctuations of their Fourier components of different frequency. Ob- 
viously such fields can be generated by illuminating the setup in Fig. 177b with a wideband 
partially coherent plane wave as described in Eqs. ()2(i4j) - ()2(i8j) . 

3. FWM's with angular noncorrelation 

According to section lVB 3l the only difference between the fully coherent FWM's and this 
special case is the lack of correlation between the tilted pulses that compose the coherent 
FWM's. Thus, to generate a FWM with angular noncorrelation one has to illuminate 
the setup by coherent light and somehow break the correlation between the tilted pulses 
propagating at different pol ar an gles. This can be done by means of the modified setup 
depicted in Fig. ED (see Ref. |l49j | for the description of corresponding quasi-monochromatic 
case). In this setup the pair of lenses L2, L3 is inserted into the path of FWM's so that in the 
focal plane of L2 the two-dimensional Fourier transform of the FWM's - the characteristic 
concentric rings - appear. In this plane we can insert a weak diffuser as to modify the 
amplitude and phase of the tilted pulses. 
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FIG. 32: Optical generation of LW's that contain plane waves that subtend non-paraxial angles 
with respect to the optical axis. 

4. FWM's with frequency and angular noncorrelation 

In this special case one has to illuminate the setup in Fig |^ with the white-light source 
as described in Eqs. - (OHH|l . 

E. Conclusions. Optical generation of general LW's 

As to conclude this chapter we note that except for the ultra-wide bandwidths, required 
to generate the wave fields described in Chapter IIVI the concept of propagation-invariance 
is well realizable in optical domain - we have shown, that the angular dispersion of various 
Bessel beam generators can be used to transform the support of the angular spectrum of 
plane waves so that to approximate that of the FWM's in some limited near-axis volume. 
The choice of the combination of the elements may be a problem in some cases, however, 
the chances are good for finding satisfactory combination. Consequently, one can launch the 
wave fields with the characteristic central peak that propagates over reasonable distances. 

Note, that the LW's containing plane waves that propagate at nonparaxial relative to z- 
axis can in principle be generated within the framework of this approach. Nevertheless, this 
requires very non- conventional optical elements like conical mirrors and diffraction gratings 
as sketched in Fig. 1221 Also, in this case the propagation distance of the generated LW's is 
quite short as can be seen from Eq. ()262p . 
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F. On the physical nature of propagation-invariance of pulsed wave fields 



As to conclude the theoretical part of this review we inter pret s ome of the properties of the 



LW's in the context of classical diffraction theory (see Ref. 113l | for a relevant discussion). 

The first note has to be made on the definition of the propagation length of LW's. Namely, 
it has been claimed in several publications (see Ref. j^i^l for example) that the LW's propagate 
over extended distances as compared to that, defined by the Rayleigh range Zji, the well 
known estimate for the scale length of the fal loff i n intensity behind of a Gaussian aperture 
in the diffraction theory, defined as (see Ref. 16(| for example) 

Zn = ^, (269) 

Wq being the minimum spot size (radius) of the beam. Indeed, for the apertured FWM's 
the radial diameter of its central spot is approximately ()75|1 

^ 4.81 _ 

kq sin Up {kq) 

(again, 7 = 1, /5 = 40^) and if we consider a planar source with diameter D = 1cm, this 
spot travels ()262|) 

h (k^in) = ^. f., = 0.9m, (270) 
2 tan [Of (fcmm)J 

The Rayleigh range ()269|) of the Gaussian pulse with radius Wq = d is Zji ^ 8.8cm ^ 
^1 (^min)- However, in our opinion such estimates are misleading and should not be used 
without the following additional details. 

The effect of diffraction is typically manifested when an obstacle is placed in the path of a 
light field. For (pulsed) beams the definition ()2(j9j) determines minimal spread corresponding 
to the waist radius Wq. In other words, this parameter determines the minimum radius of a 
circular aperture that can be placed in the path of the beam without significantly distorting 
its behaviour behind the aperture. Now, as to compare the two estimates we have to ask, 
what is the waist radius of an apertured FWM's? 

It appears, that the situation is similar to that with the monochromatic Bessel beams 
for which the energy content in every transverse lobe is appr oximately constant and equal 
to the energy content in the central maximum (see Ref. for example). In the case of 
apertured FWM's the amplitude of the branches of the characteristic X-shape fall off ap- 
proximately as 1/p. Thus, the integrated energy density as the function of radial coordinate 
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FIG. 33: The comparison of the radial field distribution of a FWM (left pane) with the field that 
is integrated over the polar angle (right pane). 



FIG. 34: The comparison of the focusing of (a) FWM's and (b) Gaussian pulses (see text). 

p is approximately constant, as the integration over the polar angle adds the factor p to 
the amplitude of the wave field (see Fig. IH^ . Thus, only a minor part of the energy of the 
apertured FWM's is contained in its central lobe and if one has to compare the propagation 
length of the apertured FWM's with the Rayleigh range of the Gaussian pulses, one should 
take Wq = D instead oiWQ = d in (IM) . 

As the matter of fact, such comparison of the Gaussian pulse and apertured FWM's is not 
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appropriate as focusing of the two wave fields are of qualitatively different physical nature. 
Namely, the Gaussian pulses are composed of monochromatic components with curved phase 
fronts, the phase fronts of the monochromatic components of the FWM's are conical. In 
lights of this difference, one can say that the FWM's are never focused in the conventional 
sense of the term and the term, focus wave mode, is rather misleading. 

In focusing Gaussian pulses most of its energy content can be concentrated into a single 
spot for a time moment. A cunning mental picture of the energetic propertied of focusing 
in FWM's can be acquired if we suppose that we have an ideal (apertured) FWM generator 
and suppose that we illuminate it with a plane wave pulse. Then central peak of each 
monochromatic Bessel beam component of the generated FWM's has the amplitude equal 
to the amplitude of the corresponding monochromatic plane wave component of the initial 
plane wave pulse. Indeed, the effect of the FWM generator to this monochromatic plane 
wave component is that its initial energy in the k-space (a finite width spot on the kz-Bjds) 
is smeared over the finite width toroid in the /c-space. In real space each point in this toroid 
will compose an apertured plane wave in the Bessel beam and the net on-axis amplitude 
of the Bessel beam is the integrated amplitude over the toroid in /c-space, i.e., again the 
amplitude of the initial monochromatic plane wave component. Now, for the superposition 
of the Bessel beams we can say, that the only place where the constructive interference 
occurs is the central spot of the apertured FWM. Thus, we can say, that given the initial 
plane wave pulse with the amplitude A and aperture D, the amplitude of the central spot of 
an apertured FWM behind an ideal FWM generator is that of the initial plane wave pulse 
A, the rest of the energy is in the sidelobes of the generated FWM's. This consequence is a 
good illustration of the qualitative difference between the FWM' and Gaussian pulse - due 
to the curved phase fronts of its monochromatic components the latter can indeed effectively 
transfer most of its energy to a single spot for a time moment. 

VII. THE EXPERIMENTS 

In our experiments we realized the optical setups for two LW's - the apertured Bessel- 

X pulses and apertured FWM's and used an interferometric cross-correlation method with 
time-integrated intensity recording to study the generated wave fields. The material of this 
chapter is published in Refs. I and V. 
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A. FWM's in interferometric experiments 



A straightforward method for recording the comphcated field shape of a coherent FWM 
or Bessel-X pulses would use a CCD camera with a gate in front of it, which should possess a 
temporal resolution and a variable firing delay both in subfemtosecond range. As such a gate 
is not realizable, any wor k able idea of experiment has to resort to a field cross- correlation 



technique (see Refs. 152| |- 157l | and references therein). 

In our experiments the wave field under investigation F (r, t) interfere with a reference 
wave Vp: 

V^{v,t) = F{v,t) + Vp{v,t). (271) 
For the reference wave we can write 

Vp{r,t)= dks{k)vp{k)exp[ik{z - c{t + At))], (272) 
Jo 

where s {k) is the (generally stochastic) frequency spectrum of the light source, vp {k) is 
the spectral phase shift introduced by the optics in the reference arm of the interferometer, 
|t;p (fc)| = 1 and At denote the variable time delay between the signal and reference wave 
fields. For the wave field under investigation we have 

/>oo 

F{r,t)= dks{k) Vp (k) Jo [kx sin Op (k)] exp[ik{z cos Op (k) - ct)], (273) 
Jo 

where vp{k) is again the undesirable spectral phase shift from the setup, |t;ir(A;)| = 1 
(we have assumed here that the FWM generator transform the input light so that for every 
spectral component s (k) the amplitude of the central spot of the corresponding Bessel beam 
is also s{k)). The averaged intensity of the resulting wave field can be expressed as 

(V^Vs) = {VpVp) + {F*F) + 2 Re {F*Vp) (274) 

(here the exact meaning of the angle brackets depends on the statistical properties of the 
light source of the experiment). 

The quantities (VpVp) and {F*F) denote time-independent intensity of the wave field. 
Specifically, the first term in the sum is the uniform intensity of the plane wave pulse: 



oo 



{V^Vp) = / dkS (k) , (275) 
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where again S (k) = {s*{k)s{k)) is the spectral density. The second term is the time- 
averaged intensity of F: 

POO 

{F*F) = / dkS (k) Jl [kx sin Op {k)] . (276) 

In principle the two components can be eliminated from the results by recording them 
separately and by numerically subtracting them from the interferograms. 
From Eqs. ()272|) and ()273|) we can write 

2Re {F*Vp) = 2Re(^J^ dkis* (ki) v*p (ki) 

X Jq (fcipsin 9f (fc)) exp [—iki {z cos 9p {k) — ct)] 

X dk2s{k2)vp{k2)exp[ik2{z-c{t-At))]Y (277) 

In our discussion we concentrated on two limiting special cases of the classification in Sec. lVBl 
- the fully coherent LW's fsection fV B 1|) and the LW's with frequency noncorrelation (section 
IV B 2|) . In both cases the averaging in (j277|) yields 

POO 

{F*Vp) = dkS (k) v*p (k) vp (k) 
Jo 

X Jo {kxsmep (A;)) exp [ikz {cos Op {k) - 1) + ikcAt] . (278) 
Here for the partially coherent field we have used Eq. (|267|) : 

{s*{k,)s{k2))=S{k,)5{k,-k2), (279) 

for coherent fields the 5 {k2 — ki) appears as the time-averaging over the term 
exp [ict {k2 — ki)]. Equation ()278|) can be given the form 



(F* Vp) = exp [i2l3-iz\ / dk S (k) v*p (ki) vp (k) 

Jo 

X Jo{kpsmeF{k))exp{ik[z{l--f) - cAt]} . (280) 

The above mathematical description yields identical results for the coherent and par- 
tially coherent fields, i.e., the results of such experiments generally do not depend on the 
correlations between the Fourier components of different frequency of the wave fields - it is 
well known that in any interferometric experiment the phase information of wave fields is 
necessarily lost. In other words, the results of the experiments do not depend on whether 
we use the transform-limited femtosecond pulses or a source of a stationary white noise. 
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The latter consequence is of great practical significance. In our overview in Chapter IIVI 
we used a spectrum that corresponds to a 3/s laser pulse and showed that the corresponding 
FWM's and Bessel-X pulse had a good spatial localization (see Fig. IT^andfTTf). However, 
computer simulations, or even simple geometrical estimations show that if the autocorrela- 
tion time of the source field r exceeds ~ 10 femtoseconds, the characteristic X branching 
occurs too far from the axis z and in this narrow-band limit the resulting wave field would 
be nothing but a trivial interference of quasi-monochromatic plane waves. Thus, the band- 
width of the light source is a very challenging part of the setup. In what follows we add to 
the reputation of incoherent sources as being the poor mans femtosecond source and confine 
ourselves to the special case of frequency non-correlating fields. 

B. Experiment on optical Bessel-X pulses 



Our setup for the interferometric experiments on optical Bessel-X pulses is depicted in 
Fig. IH^b - as compared to the Bessel-X pulse generator in Fig. 1221 the pinhole is made in 
the centre of the annular ring mask as to form the reference field (plane wave pulse) behind 
the lens L. For the mathematical description of the situation we have to choose /5 = in 
Eq. ()280j) . i.e., d (k) = 9q = arccos7 = arccosc/f . Setting also At = we get 



where S (k) is the spectral density of the light source (note that the annular ring mask 
has the uniform spectral response function as explained in the end of the Sec. IVl B|) . In 
Eq. ()281|1 Zm denotes the distance along the optical axis of the setup. To understand the 
significance of this parameter we have to remember the superluminal group velocity of the 
Bessel-X pulse - during the flight the latter catch-up with the reference plane wave pulse. 
In this context the coordinate Zm is the distance of the recording device from this catch-up 
point (see Fig. HEb). 

If we define in the general expression for the mutual coherence function of the wave fields 
in Eq. ()21fij) the origin of the z axis as being in the point z = ct so that Az can be replaced 



1. Setup 




X Jo {kp sin ^o) exp [ikz^ (1 — 7)] 



(281) 
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FIG. 35: (a) Intensity profile of a computer-simulated Bessel-X pulse flying in space, shown as 
surfaces on which the field intensity is equal to a fraction 0.13 (l/e^) of its maximum value in 
the central point. The field intensity outside the central bright spot has been multiplyed by the 
radial distance in order to reveal the weak off-axis side-lobes. Inset: amplitude distribution in the 
plane shown as intersecting the pulse. The plots have been computed for 3-fsec near-Gaussian- 
spectrum source pulse with carrier wavelenght A = 0.6^m and the angle 9o = 0.223 deg. For these 
parameters the dimensions of the plot xyz box are 20 times 20 times /um; (b) Optical scheme of 
the experiment. Mutual instantaneous placement of the Bessel-X pulse and the plane wave pulse 
is shown for three recording positions (two of which labeled in accordance with Fig. I36|) . The 
ovals indicate toroid-like correlation volumes where co-propagating Bessel-X and plane wave pulses 
interfere at different propagation distances along the z-axis. L's, lenses; M, mask with Durnin's 
annular slit and an additional central pinhole for creating the plane wave; PH, cooled pinhole 10/xm 
in diameter to assure the transversal coherence of the light from the source V. In case source V 
generates a non-transform-limited pulses or a white cw noise, the bright shapes depict propagation 

of the correlation functions instead of the pulses. 

Ill 



by the distance from the pulse centre zi^ and set r = 0, r_Li = r^, rj_2 = the result reads 
rF(r±,0,ZA) =exp[i2/372A] / dfc |V (A;)|^ Jq (A;psin^F (A;)) exp [2fc7ZA] • (282) 
If we also set /? = we get 

/•oo 

Tbx (r±, 0, za) = / dA; |V (A;)|^ Jo (^psin^o) exp [iA;72;A] . (283) 

Comparing the Eqs. (j2HIl) and (j2HSl) we see that if |V {k)f = S (k) and v*p (ki) vp (k) = 1, 
we have 

{F*Vp) = Tbx 0, z^^—^] , (284) 



7 

so that 

ZA = zj—^. (285) 

7 

The interpretation of the small factor (1 — 7) /7 in ()285|1 is that the setup serves as a "z 
axis microscope" for recording the mutual coherence function 1)2811) along the 2;-axis which 
scales the micrometer-range z dependence of the field into a centimeter range. 

If we compare the expression ()283p with the hypothetically measurable field distribution 
given as the real part of the Bessel-X pulse in Eq. ()135p . 

/•oo 

^ Bx (p, z,t) = dk s (k) Jo [kp sin 6^0] exp [ik {z-f - ct)] (286) 
Jo 

we conclude that the experiment reveals the whole spatiotemporal structure of the Bessel-X 
field. The natural price we have to pay for resorting to the correlation measurements is 
replacing the spectrum s (k) with its autocorrelation, which is a minor issue in the case 
of transform-limited source pulses. Nevertheless, we cannot claim, that we actually detect 
the field under investigation (see also Ref. 0] for a relevant discussion). Indeed, as the 
absolute phases of the plane wave components are inevitably lost in any linear interferometric 
experiments only the spatial amplitude distribution of the wave field can be detected. 

In the reasons described above we took advantage of the insensitivity of Eq. ()28ip to 
the source field phase and used a white light noise from a superhigh pressure Xe-arc lamp 
instead of a laser as the field source to achieve the ~ 3 femtosecond correlation time in our 
experiment {V in Fig. IHKb). Thus we implemented the third special case of Sec. IV Bl the 
Bessel-X field with frequency noncorrelation. 
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FIG. 36: (a) Samples of the experimental recordings and processing of the intensity distributions 
measured at the positions Pos.A and Pos.B along the z-axis as shown in Fig. 135b . Left column, 
total interference pattern of the cross-correlated fields; middle column, lateral intensity distribution 
of the Bessel-X field alone. In the right column gray (about 50 per cent) corresponds to zero level 
and dark to negative values; (b) Comparison of the result of the experiment (right panel) with a 
computer-simulated Bessel-X pulse field. 

2. Results of the experiment 

The recordings at 70 points on the z axis (from behind the L3 lens up to a point a 
few centimeters beyond the origin) with a 0.5-cm step were performed with a cooled CCD- 
camera EDC-IOOOTE, which has 2.64 x 2.64 mm working area containing 192 x 165 pixels, 
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and processed by a PC as follows. 

First, the subtraction of the Bessel-X field intensity was performed (see Fig.l^Hk). whereas 
the same procedure with the plane wave field intensity, due to its practically even distribu- 
tion, was found to be unnecessary. In order to reduce noise and the dimensionality of the 
data array, an averaging over the polar angle in every recording was carried out by taking 
advantage of the axial symmetry of the field. Thus we got a ID array containing up to 
hundred significant elements from every 192 x 165 matrix recorded. Seventy such arrays 
formed a matrix, which, having in mind the known symmetry of the real part of Eq. ()28H) . 
was mirrored in the lateral and the axial plane. The result is compared in Fig. IHHb with 
the Bessel-X field distribution in an axial plane, computed from Eq. (j286|) for a model 
spectrum. The latter was taken as a convex curve covering the whole visible region from 
blue to near infrared (up to 0.9 fim) in order to simulate the effective light spectrum in 
the experiment, which is a product of the Xe-arc spectrum and the sensitivity curve of the 
camera. The central (carrier) frequency was chosen corresponding to wavelength 0.6 fim 
which had been determined from the fringe spacing of an autocorrelation pattern recorded 
for the light source with the same CCD camera. The left- and right-hand tails of the central 
X-like structure are more conspicuous in the experimental pattern due to unevenness of the 
real Xe-arc spectrum. The different scaling of the horizontal axes of the two panels is in 
accordance with the z axis "magnification" factor of the experimental setup. 

Observing the obvious agreement between theoretical and experimental patterns, we ar- 
rive at a conclusion that we have really recorded the characteristic spatiotemporal profile of 
an optical realization of the nonspreading axisymmetric Bessel-X field. 

C. Experiment on optical FWM's 

1. 3D FWM's and 2D FWM's, the mathematical description of the experiment 

To prove the feasibility of the approach for optical generation of apertured FWM's de- 
scribed in Chapter |VT1 one has to implement the setup depicted in Fig. |2Zb and show that 
the generated wave field indeed behaves as a FWM in interferometric experiments. However, 
without the loss of generality the task can be simplified as follows. 

In section IIII C II we demonstrated that all the defining properties of the FWM's - its 
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propagation-invariance, the characteristic field distribution and pre-determined group veloc- 
ity - can be studied in terms of the specific pair of interfering tilted pulses, the 2D FWM's 
in Eq. 



F2D{x,y,z,t)= / dkB{k) 
Jo 

X COS [kx sin 9f (k)] exp [ik {z cos 6f (k) — ct)] (287) 

(here we have set = 0). In other words, we have shown that the peculiar propagation of 
FWM's is assured exclusively by the specific coupling between the wave number and the 
direction of propagation of plane wave components of the FWM's as defined by the function 
Opik) in Eq. 



In 56| we experimentally demonstrated the feasibility of 2D FWM's in Eq. (128 7|) as the 
physical concept which is much more transparent in this sort of experiments. However, it 
also appeared, that the fabrication and polishing of a high-quality, large aperture concave 
conical surfaces is still a complicated task. 

The mutual coherence function for the 2D FWM's can be easily deduced from the one 
for the (3D) FWM's ()21fj|l . In complete analogy with the discussion in section IlII C II we 
replace the Bessel function Jq () by cos (), choose appropriate coordinates and get 

poo 

r2Z) (a^i, a;2, A2;, 7Az — cr) = exp [— Z/37A2;] / dk \V [k)^ 

Jo 

X COS {kxi sin Op (k)) cos {kx2 sin Op (k)) exp [ik {•jAz — cr)] . (288) 

Also, the mathematical description of the experiments with 2D FWM's is analogous to that 
of the three-dimensional one. In the spatial intensity distribution of the interferometric 
experiment 

{V*Vj:) = {V^Vp) + {F;^F2d) + 2 Re {F*^Vp) , (289) 
we have for the intensity of the 2D FWM 

POO 

{F2DF2D) = / dkS (k) cos^ {kx sin Og (k)) (290) 
Jo 

where the angular function 9g {k) [see Eq. ()242|) ] is determined by the specific setup and 
S {k) is the spectral density of the light source (see notes in the end of the Sec. IVI B|) . For 
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the third term instead of Eq. ()28Up , we have 

poo 

{F;dVp) = exp [2p^z] / dA; S (k) v*p (ki) vp (k) 

Jo 

X COS {kp sin edk)) exp {ik[z{l--f) - cAt]} . (291) 
In P7T|1 for FWM's we have to set 7 = 1, so that 

POD 

{F*jjVp) = exp [2pz] / dkS (k) v} (ki) vp (k) 
Jo 

X COS {kp sin 9g (k)) exp [—ikcAt] . (292) 



In analogy with the case of Bessel-X fields we can define in ()288|) Az = 0, r_|_2 = 0, so 
that the mutual coherence functions of the 2D FWM reads 

POO 

TsD (x, 0, 0, -cr) = / dA; |V(A;)|^cos(A;xsin^G(A;))exp[-iA;cr] . (293) 
^0 



Thus, if we record the interference pattern at 2; = 0, the comparison of Eqs. ()293|) and 
yields 

{F2DVp)^t = T2D (x, 0, 0, -cAt) , (294) 

and we can conclude that the mutual coherence function of the 2D FWM's can be studied 
by recording the intensity of the interference picture as the function of the delay At between 
the signal and reference fields. 

Note, that the equation ()292|) can be given the form 



POD 

{f;^Vp)= / dksik) 

Jo 

X COS {kx sin 9g (k)) exp [ikz cos Oq (k) — ike (to — At)] , (295) 

where to = z/c and the constant has the interpretation of being the time that a wave field 
propagating at group velocity c travels the distance z to the plane of measurement. The 
integral expression ()295p is very similar to the one describing the field of the 2D FWM's 
(128 7|) . the only difference being that the theoretical angular function Op (k) is replaced 
by the Ocik) and frequency spectrum is replaced by the power spectrum S {k) in ()295j) . 
Again, as the absolute phases of the plane wave components are inevitably lost in any 
linear interferometric experiments we can detect only the amplitude distribution of the wave 
field. However, the interferograms do carry information about the defining, most essential 
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characteristic of the FWM's - their support of the angular spectrum of plane waves. Indeed, 
the general structure of the interference patterns in Eq. ()295|) is primarily determined by 
the angular function 6g (k), and it resembles the corresponding transform-limited wave field 
only if 6g (k) = Op {k) in good approximation over the entire bandwidth of the field. 



2. Setup 

The conversion of the FWM generator in Fig. HTh to the two-dimensional case is straight- 
forward - we just replace the axicon and circular diffraction by their one-dimensional coun- 
terparts - prisms (wedges) and a diffraction grating. The initial field on the elements is an 
interfering pair of pulsed plane waves. 

The setup of our experiment is depicted in Fig. EZI The main part of it is the (2D) FWM 
generator that consist of the mirrors M7 and M8, of the two wedges Wl and W2 and of a 
blazed diffraction grating G.(see grayed area in Fig.EZj). The FWM generator is placed into 
an arm of an interferometer as will be explained below. 

In our experiment we implemented a 2D FWM with following parameters: f3 = AOrad/m, 
j = l {vg = c) giving Of (ko) ^ 0.23° if ko = 7.8 x lO^rad/m (Aq ^ SOOnm) (see Fig.EHt for 
the spectral density of the light source and Fig. IHHb for the support of the angular spectrum 
of plane waves of the specified 2D FWM). 

The FWM generator has three free parameters: a - the apex angle of the wedges, ^ 
the angle that the initial pulsed plane waves subtend with the optical axis of the setup and 
d - the groove spacing of the diffraction grating. As to find the values for the parameters 
that give the best fit between the generated support of the angular spectrum of plane waves 
and the support of angular spectrum of plane waves of the theoretical FWM's in Eq. (jMj) 

9f (k) = arccos ( (296) 



k 

(7 = 1) we combine the ()296p with the Eq. ()243p and write the following system of equations: 
arccos / li^^nM J = ^ + a {1 - n (fc„)) + 9o,m = l,2,3. (297) 

We specify the three wave numbers as ki = 7.4 x 10^^ (A = 849nm), ^2 = 10.6 x 

(A = 593?T,m), ks = 1.6 X 10^^ (A = 381nm) and assumed that the wedges are made of the 

optical glass BK 7 for which the refractive index n (k) is known to an accuracy better than 
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FIG. 37: Experimental setup for generating 2D FWM's and recording its interference with plane 
wave pulses. The FWM generator can be seen in grayed area; M's, mirrors; L's, lenses; BS's, beam 
splitters; W's, wedges; G, diffractional grating; AL, Xe-arc lamp; PH, pinhole; GP's, compensating 
glass plates. 

10"^ The system pJ7j) then yields 

a = 1.2044 X lO'Vad 

d = 3.7494 X 10~^m (298) 
60 = 9.4683 X IQ-^rad . 

As inserted into Eq. ()242p . the maximum deviation 69 {k) = Op {k)—6c {k) for the wavelength 
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FIG. 38: (a) The power spectrum of the hght used in our setup; (b) The angular spectrum of the 
plane waves generated in the setup (solid black line) as compared to the theoretical one (dotted 
cyan line); (c) The deviation of generated support, here ki = 7.4 x 10^^ (A = 849nm), /c2 = 
10.6 X lO'^^ (A = 593nm), = 1.6 x lO"^^ (A = 381nm). 

dispersion of the cone angle in the selected spectral range is as small as 5 x 10~^deg, i.e., 
< 0.2%. A comparison of the corresponding supports of angular spectrum and the exact 
form of deviation 69 {k) are depicted in Figs. IHHb and IHHh (see also Fig. 1^ . As for rough 
estimation of the spread of the pulse due to the 56 {k) one can estimate the corresponding 
maximum group velocity dispersion Av^ and compare this with the mean wavelength of the 
pulse. The numerical simulations show that the approximation is indeed good enough for 
propagating the central peak of the FWM's over several meters. 
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The FWM generator has been placed in what is basically a specially designed, modi- 
fied Mach-Zehnder interferometer. The interferometer consists of two beamsplitters and of 
identical broadband non-dispersive mirrors. The input field from the light source is split by 
the beamsplitter BSl into the fields that travel through the two arms of the interferometer, 
the one with the FWM generator and the arm for the reference beam. The mirrors M5 
and M6 form a delay line, they were translated by the Burleigh Inchworm linear step mo- 
tor, the Ifim translation step of which was reduced to Qbnm by a transmission mechanism. 
The mirror M7 was continuously translatable as to correct for the time-shift between the 
two tilted pulses. The wedges Wl and W2 were transversally translatable as to balance 
the material dispersion they introduce to the plane wave pulses (see text below). We used 
Kodak Megaplus 1.6i CCD camera with the 1534 x 1024 matrix resolution and 10 bit pixel 
depth. The linear dimensions of the matrix are 13.8mm(H) x9.2mm(V), the pixel size is 
9/xm X 9/xm. 

Again, in our experiment we used the filtered light form a superhigh pressure Xe-arc lamp, 
giving ^ 5 fsec correlation time for the input field [see Fig. IHHk for the power spectrum of 
the light]. To ensure good transversal coherence over the clear aperture of the setup the 
required maximum diameter 15/im of the pinhole and focal length 2m of the collimating 
Fourier lens LI was estimated from the van Cittert-Zernike theorem jlGS^ for the mean 
wavelength of the light Aq = SOOram. As the result of filtering, the total power of the signal 
on the ~ 1.5cm^ CCD chip was very low, approximately 0.03/iiy. 

Due to the short coherence time of the source field, the experiment is highly sensitive to 
the phase distortions (spectral phase shift) introduced by the dispersive optical elements of 
the system - the beamsplitters and the FWM generator. In the FWM generator there is 
three possible sources of undesirable dispersion: (1) the propagation in the glass substrate 
of the diffraction grating, (2) the propagation between the grating and the axicon where 
the support of angular spectrum of the wave field is not appropriate for the free space 
propagation, i.e., it does not obey the Eq. (jH^ and (3) the propagation in the wedges. The 
beamsplitters in our setup are identical and if we set them perpendicularly and orient the 
coated sides so that each beam passes the glass substrate of the beamsplitter twice, the 
arms of the interferometer remain balanced. The infiuence of the propagation between the 
elements can be made neglible by placing them close to each other. The character of the 
undesirable dispersion in the wedges can be estimated from the following considerations. The 
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entrance wave field on the wedges is the transform-limited Bessel-X pulse, so, the on-axis part 
of the pulse is also transform-limited and should pass through the wedges unchanged, i.e., 
without any additional spectral phase shift. Consequently, the wedges should be produced 
and aligned so that their thickness is zero on the optical axis. As the apex angle of the wedges 
is very small in our setup (~ 0.7°), this is not very practical approach and we consider the 
finite thickness on the axis as the source of additional spectral phase shift instead. Thus, 
the composite spectral phase shift of the FWM generator can be described as the phase 
distortion introduced by the substrate of the diffraction grating and by a glass plate of the 
material of the wedges, the thickness of which is equal to the thickness of the wedges on the 
optical axis. We balanced the arms of the interferometer by inserting material dispersion 
into the reference arm of the setup by means of two appropriate glass plates (GPl and GP2 
in Fig.EZK). 

3. Results of the experiment 

In first experiment we recorded the time-averaged interference pattern of the 2D FWM 
and the reference wave field as the function of the time delay between the two. The exper- 
iment can be mathematically modeled by varying parameter At in Eq. ()292|) . We scanned 
the time-delay at three z-axis positions, z = 0cm, z = 25cm, z = 50cm (the origin of the 
z-axis is about 30cm away from the beamsplitter BS2 in Fig. ITfk ) . In each experiment we 
recorded 300 interferograms, the time-delay step was 0.43 /sec (0.13 /im). 

In a typical interference pattern in our experiment [see Fig. IHIIk ] the sharp vertical inter- 
ference fringes in the center correspond to the second term in the interference sum fl289|) - 
this is the time-averaged, "propagation-invariant" time-averaged intensity of the 2D FWM. 
The fringes can also be interpreted as the autocorrelation function of the interfering tilted 
pulses [see Fig. IHHk for the corresponding power spectrum] . In this experiment the intensity 
of the wave field under study do not carry any important information, so we subtracted it nu- 
merically from the results in Fig. 139b. The interference fringes that are symmetrical at both 
sides of the central part correspond to the third, most important term in this sum. It can 
be seen from Fig. l39k . that due to the low signal level the recorded interferograms are quite 
noisy. To get the better signal-to-noise ratio, we averaged the data in the interferograms 
over the rows and used the resulting one-dimensional data arrays instead. 
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FIG. 39: (a) A typical interferogram in the setup as recorded by the CCD camera; (b) The 
interference pattern in the setup as the function of the delay between the signal and reference wave 
fields in three positions of the CCD camera (see text for more detailed description). 

The results of the experiment are depicted in Fig. I^Ub . We can see, that there is a good 
qualitative resemblance between the measured xAt plot of the interference pattern and the 
theoretical field distribution of the 2D FWM's in Fig. [7|d as it was predicted by the Eq. (j295p 
- one can clearly recognize the two interfering tilted pulses forming the characteristic X- 
branching, also the phase fronts in the tilted pulses can be seen. The wave field is definitely 
transform-limited, so we have managed to compensate for the spectral phase shift in the 2D 
FWM generator. 

We can also see, that the interference pattern does not show any spread over the 0.5 m 
distance, consequently, the wave field does not spread in the course of propagation. 
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FIG. 40: The interference pattern as the function of the CCD camera position (see text for detailed 
description). 

An additional detail can be found in Fig. 1^% : the tilted pulses do not extend across 
the whole picture but are cut out (see dashed lines in Fig. Also, the "edges" of the 

tilted pulses move away from the optical axis. This effect can be clearly interpreted as the 
consequence of the finite extent of the tilted pulses, as illustrated in Figs. ElK andEEl^ the 
dashed lines simply mark the borders of the volume, where the two tilted pulses intersect, 
i.e., the borders of the volume, where the 2D FWM exist [see the striped area in Fig. Elk]. 
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In the second experiment we recorded the interference pattern as the function of the prop- 
agation distance z. The experiment can be simulated by varying z coordinate in Eq. ()292p . 
We recorded 240 interferograms, the step of the CCD camera position was 3.1mm. The 
numerical simulation of the experiment and the results of the experiment are depicted in 
Figs. liUb and liUh respectively. 

The experiment can be easily interpreted - the position-invariant envelope of the inter- 
ference pattern is the consequence of the fact, that the group velocities of the propagation- 
invariant 2D FWM and the reference field are equal, c, so that the overlapping volume of 
the two fields do not change in the course of propagation (see Fig. HHk). The z dependent 
finer structure of the interferograms is the consequence of the fact, that the phase velocities 
of the plane wave pulse and 2D FWM are not equal, i.e., we have also 7^ for the phase 
and group velocities of the 2D FWM. The result of the experiment in Fig. HUb show good 
qualitative agreement with the theory. 

We can also determine the parameter /? from our experiment - the exponent multiplier 
in Eq. ()292j) reads exp [i2/3z], thus /5 = tt/zq where zq is the period of the variations along 
the 2;-axis. From the result in Fig. HHk we estimated zq ~ 7.5cm, so that /3 ~ 42 rati/m, 
which result is in good agreement with the theory. 

Thus, we have shown, that the generated wave field has all the characteristic properties 
described in previous theoretical sections and the validity of the general idea has been given 
an experimental proof. 



VIII. SELF-IMAGING OF PULSED WAVE 
FIELDS 



Self-imaging, also known as Talbot effect is, in its original sense, a well-known phe- 
nomenon in classical wave optics where certain wave fields reproduce their transvers al in- 
te nsity distribution at periodic spatial intervals in the course of propagation (see Refs. \l6^ 



177| . and references therein). The effect has been studied extensively by means of Fresnel 



diffraction theory and the angular spectrum representation of scalar wave fields. As a result, 
the general, physically transparent conditions have been formulated the transversal intensity 
distribution of a wave field have to obey to be self-imaging. It is also well known that the 
monochromatic propagation-invariant wave fields (Bessel beams) constitute a special class 
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of self-imaging wave fields, the mathematical description of the two being much the same. 
In recent years the effect has attracted a renewed interest, as the concept has been 



generalized in to the domain of pulsed wave fields by several authors (see Refs. 



176 



18G 



181 



172 
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, and references therein) . The phenomenon has been discus sed in the context o f 



fiber optics jlT^ and also as a property of spatial, wideband wave fields 
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l3 
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In what follows we show, that the discrete superpositions of FWM's over the parameter 



/9 



(299) 



can be used as the self-imaging "pixels" of a spatiotemporally self-imaging three-dimensional 
spatial image - the wave fields of this type can reproduce spatial separated copies of its initial 
three-dimensional intensity distribution at specific time intervals. The results in this chapter 



are published in Ref. 



55|. 



(I ; is important to note, that our discussion is closely related to those in Refs. [ 



and 18 1| - those publications consider essentially the same problem, however the analysis 
is different in each occasion) 



A. Monochromatic self-imaging 



The se lf-irn aging condition for the monochromatic wave field has been defined as (see, 
e.g., Ref. |l3) 

^{p,z,^,t) = ^>{p,z + d,^,t). (300) 

With the condition ()300|) and the Whittaker type angular spectrum of plane waves of 
monochromatic wave fields (see Eq. (fTTj) 



\1/ (p, z, ip, t; k) = exp [±m(y9] / d9 {k, 9) 
n=o 



X Jn [kp sin 0) exp [ik {z cos 6 — ct)] 



(301) 



one can eas ily d educe the condition for self-imaging for the monochromatic scalar wave fields 
that reads |l69l | 



kd cos 9 = + 27rg, 



(302) 
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where q is an integer and ip is an arbitrary phase factor. The relation ()302p imphes, that a 
monochromatic wave field ()15|) periodically reconstructs its initial transversal field distribu- 
tion if only its angular spectrum of plane waves is sampled so that the condition 

%!) + 27rg 



d 



(303) 



where is again the z component of the wave vector, is satisfied for every plane wave 
component of the wave field. Applying the condition ()303|) and including only the axially 
symmetric terms in the summation, we get the following expression for the general cylindri- 
cally symmetric, monochromatic, self-imaging wave field: 



\1/ (p, t; k) = exp [— icut] a, Jq 
where we have denoted 




exp 



'.27Tq ' 
I — —z 



A^ik, arccos 



27rg 
kd 



(304) 



(305) 



and '?/' = is assumed. Thus, we have a discrete superposition of Bessel beams the cone 
angles are specified by Eq. ()303p . The physical content of Eqs. ()303p and ()304|) can be 
summarized by saying that the monochromatic self-imaging is essentially the phenomenon 
where the wave field is a discrete superposition of wave fields with different phase velocities 
so that the total field, being the superposition of the composite fields, depends periodically 
on the distance. 

The on-axis superposition in 



\1/ (p, z, t; k) = exp [— icut] exp 



■.27rg ■ 



(306) 



can be recognized as the Fourier series representation of a periodic function along the optical 
axis. However, in ()306|) the condition ()302j) together with the causality requirement k^ > 
imply that 



< g < 



kd 
2^' 



(307) 



An example of the on axis intensity distribution of a superposition of Bessel beams in ()3()6|l 
is depicted in Fig. 



126 



FIG. 41: The snapshots of the temporal evolution of the spatial intensity distributions of the 
pulsed self-imaging wave field, the three-dimensional self-imaging "pixel" , at 26 femtosecond time 
intervals. Due to the axial symmetry, the distribution is shown in one meridional (say, xz) plane. 



B. Self-imaging of pulsed wave fields 



From the material in Chapter IIIII we know that the phase velocities of the FWM's 
with different parameter f3 along the optical axis are generally different and we can write 

= cos^H^o,/?) ^ 7 (ko - 2/5) • 



Thus, the superposition of FWM's over the parameter /3 in 

(p, z, ^,t)=J2 C.^i^ i^, P. t; P,) (309) 
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is taken over a set of overlapping, non-spreading optical pulses that (1) are transversally 
localized, i.e., their transversal intensity distribution have a single narrow intense peak, (2) 
have equal carrier frequency, (3) propagate at equal group velocities, but (4) have different 
phase velocities. In complete analogy with the monochromatic self-imaging we could suggest 
that if the component pulses satisfy certain conditions, the single narrow peak of the wave 
field could periodically vanish and reconstruct its initial (localized) transversal intensity 
distribution. As a result we could get a wave field the temporal evolution of which can be 
perceived as a spatial arrow of sequentially visible light spots (see Fig. HTj) . In this sense, 
such superposition could be considered self-imaging and, due to its spatial localization, it 
could be used as a pulsed "self-imaging pixel" of a transversal or even spatial image. 
Consider a discrete superposition of a set of the tilted pulses in Eq. (jHH|) 



jn{SI) 



dA; A (k) exp [ik {x sin Op (k) + z cos 6f {k) — ct)] 



(310) 



the function 6p{k,(3q) being determined by the Eq. (jMj) . The expression can be given a 
readily interpretable form if we approximate the x and z components of the wave vector by 



k^ {k, p) = k^ {ko, p) + ^k^ {k, p) 

ak 



{k - ko 



(311) 



k=ko 



k,{k,f3) = k,{ko,(3)+ -^k,{k,f3) 

ak 

= k^ {ko,p) + 'j{k-ko) , 



{k - ko) 



k=ko 



(312) 



where we have used Eq. (jH^ . Substitution of relations ()31ip and ()312|) in Eq. ()310|) yields 

r(^^) (x, z,t)^Y^ aqLq (x, z-i - ct) 
q 

X exp [iko (x sin Op {ko, Pq) + z cos 9f {ko, Pq) — ct)] , (313) 

where exp [iko {x sin 6p {ko, Pg) + z cos Op {ko, Pq) — ct)] is the carrier wavelength plane wave 
component of the tilted pulse ()310p and 



Lq {x, Z7 - ct) 



dkA {k + ko) 



■ko 



X exp 



ik I X '^^x {k,Pq) 



+ — ct 



ko 



(314) 
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is an approximation to the non-spreading traveling envelope of the pulse. In near axis volume 
we can write 

T^^^) (x, t) = L (x, ^7 - ct) ^ flg (315) 

X exp [iko (x sin Op {ko, [3q) + z cos Op {k^, (3q) — ct)] . 

The Eq. fl315j) is essentially a product of a propagating pulse L(x, 27 — ct) and of a term, that 
is a mathematical equivalent of the superposition of the monochromatic carrier-wavelength 
plane waves that angles 9p {ko, (3g) to the optical axis. According to our general 

idea, the latter term should be self-imaging in the conventional, monochromatic sense of the 
term in Eq. ()304|) - in this case the product ()315|) behaves as a pulse that vanishes and 
reconstructs itself periodically. 

In complete analogy with the Eq. ()3()2|1 the condition for such behavior can be written as 

k = dcosOp {ko^ (3q) = tp + 2Tiq. (316) 

and with the Eq. (jHSJ we can write 

^ = 7^0 - 27/? , (317) 
so that we get a discrete set of constants {3 for the superposition ()310|) : 

"-t-^ (31«) 
(see Fig. 021 for an example). The direction of propagation of the carrier wavelength plane 
wave component of the tilted pulse ()310p can be found by combining the Eqs. (jH^ and ()318p : 
we can write 

Of {k, Pg) = arccos + '^q^ , (319) 

so that 

Of {ko, Pq) = arccos ( ^q] ■ (320) 



kod ' 

Note, that again, in ()320|) the angle Qp {k,Pq) have to satisfy the condition 

knd 

< g < ^ . (321) 
With those conditions we can write the cylindrical superposition of the tilted pulses in 

Eq. dnnD 

/■27r 

vl/(5^) (p, z,t)= / dip T(^^) (x, z, t; ^) (322) 
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FIG. 42: An example of the set of supports of the angular spectrums of plane waves of a self-imaging 
superposition of the tilted pulses (see text). 



FIG. 43: (a) The Fourier spectrum and (b) the spatial amplitude of a train of sinusoidal waves 



as 



[p, z,t) = ag exp [-i2-f(3gz] / dkA (k) 



(323) 



X Jn 



kpJl- 



k 



exp [ik (^7 — ci)] . 



and this is the discrete superposition of FWM's we were looking for in Eq. ()309p . 

The on-axis longitudinal shape of the superposition in Eq. ()323p can be easily evaluated for 
the most practical case of a uniform superposition of {2n + 1) tilted pulses, centered around 
some carrier spatial frequency kz {k^, (3q) (see Fig. HHj). Superposition can be expressed as 

sin [\ {2n + 1) Ak^z] 



^ z) 



A. 



sin {j2^kzZ^ 



(324) 
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where A/c^ is the interval between the spatial frequencies (see Ref. (l67l | for example). For 
this case the self- imaging distance is d = 2TT/Akz and the width of the peaks of the resulting 
function is Az ~ 27r/(2?2+l) A/c^. The result of the evaluation of Eq. (|32U|) for a superposition 
of seven tilted pulses is shown on Fig. I^Hb . 

A numerical example of the self-imaging behavior of a superposition of five FWM's in 
fl323p is depicted in Fig. |^ In this example the self-imaging distance d = 2x 10~^m, 7 = 1, 
ak = 3.8 X 10-^m 3/s) q = 27, 28, 31 with (3g being determined 
by Eq. {(3q = 7.6 x 10W/m...l.3 x lOW/m) and Op {ko, (3q) by Eq. (jH^ . 

As the second example we demonstrate the self-imaging transmission of a non-trivial 
spatial image, depicted on Fig. 144b. The image consists of eight "pixels"- the self-imaging 
superpositions of FWM's - specified in previous example. The numerical examples clearly 
show that the concept, in principle, is applicable for constructing wave fields that self-image 
three-dimensional images. Still, the experimental realization of such wave fields is not trivial. 



IX. CONCLUSIONS 



In this review we developed a physically transparent, comprehensive theory for the de- 
scription of propagation-invariance of the scalar wideband wave fields in terms of Whittaker 
type angular spectrum of plane waves. This representation was demonstrated to be very 
useful in discussions on the physical nature of the localized wave transmission phenomenon 
and propose a general idea for optical generation of the localized waves. 
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FIG. 44: A numerical example of the evolution of a self-imaging spatial image (smiling human face) 
consisting of eight self-imaging pixels. The snapshots are taken at 19 femtosecond time intervals. 
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X. NOTATIONS USED IN THIS REVIEW 



PiZ^i-p - cylindrical coordinates system 

ky, kz - the Cartesian components of the wave vector k 

/cq - carrier wave number 

k,9,(j) - spherical coordinates in /c-space; 
X = kp = k sin 9 

7, 7p - see Eq. §B 

(3,i - see Eqs. ^ and (pl 

fi = z + ct 
( = z — ct 

v^, - group velocity ()32j) and phase velocity 

Oi, 02, K, 6, p, g - parameters, see Eqs. ()115|) . ()118|) . ()155|) 

Op (k) = Op ik,p), e'^p (k) - see Eqs. (jSl, © 

Og (k) - see Eq. §^ 

kp (0) - see Eq. (j^ 

« = |(^ + fc.) 

/3=|(f-/c.) 

T = t2 — ti, Az - the time- and z coordinate difference 

Ts, o"z, o'p - see Eqs. (jHSj), ()66a|) . ()66b|) 

At, ZA, Zm - see Eqs. (jTZ^ . (j2H21), 

(k, cj) = (r, t)] - see Eq. © 

A (/ca;, - see Eq. (fT^ 

A„ (A;, 6) - see Eq. jlTj) 

(fc^, fcy), (k,) - see Eq. dSni), il 

A;^*^ {k, x) - see Eqs. (EH) 

Bnik), Bo{k) ^ B{k) - see Eq. (jUl) 
y4 (/ca;, kykz), An {k, 9), En {k) - see note after Eq. (jH3j) 

- see Eq. (El 

S (x, /3) - see Eq. 

a(fc, 6', 0) - stochastic angular spectrum of plane waves 
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E(r,t), H(r,t),A (r,t), U^^\ H^™) 

(p, 2:,</',i), ^/ {p,z,ip,t) 
T {x,y,z,t;(f)), F{x,z,t) 

r (ri,r2,ti,t2) 
W^(ri,r2,fc) 

v4(/Cl,fc2,6'i,6'2,01,02,) 

C {ki,k2,(f)i,(f)2) 
s{k) 

S{k) = \s{k)f 



see Eqs. ^ 



see Sec. Ill Al 



see Eqs. dH, (HH 
see Eqs. (jHHI), (ffl) 

mutual coherence function, Eq. p77|) 
cross-spectral density 
see Eq. angular correlation function 

see Eq. 



see Eq. (fTO 

see Eq. (jMj) . frequency spectrum of light source 
spectral density (power spectrum) of light source 
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